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ÜBER DIE AUSBREITUNG 
„EXPONENTIAL BEGRENZTER“ WELLEN 
VON ENDLICHER AMPLITUDE 


Von Hara RYFFERT 
Institut der Theoretischen Physik der Universität Poznań 


( Erhalten am 6. Oktober 1954) 


Es wird die Ausbreitung einer akustischen Welle in einer Exponentialbegrenzung 
untersucht, wobei das Problem mit Hilfe der krummlinigen Koordinaten als ein eindimen- 
sionales Problem formuliert wird. Die Wellen sind durch harmonische Schwingungen 
eines Kolbens, der sich im engen Teil der Begrenzung befindet, hervorgerufen. Den 
Ausgangspunkt bilden die hydrodynamischen Bewegungsgleichungen in Lagrangeschen 
Koordinaten, die dadurch linearisiert werden, dass alle Grössen nach Potenzen der Rela- 
tivamplitude A = 2ла/5 entwickelt werden, wo a die Amplitude der harmonischen Schwin- 
gungen des Kolbens ist und A die Wellenlänge, die der Kreisfrequenz o der Kolbenschwin- 
gungen entspricht. In erster und zweiter Näherung haben die entstehenden harmonischen 
Wellen die Kreisfrequenzen beziehungsweise 2c. Die Amplitude der Grundwelle ver- 
kleinert sich wie die Funktion e—9*/2, dagegen die harmonischen Wellen verkleinern sich 
teilweise wie die Funktion e—9*/2 und teilweise wie e—4*, 

Die Arbeit enthält Formeln zur Berechnung des Verhältnisses der Intensitäten dieser 
beiden Wellen an der Mündung der Begrenzung. Dieses Verhältnis erweist sich als propor- 
tional dem Quadrate der Relativamplitude A. Endlich wird der Schluss gezogen, dass 
ein Lautsprecher mit Exponentialbegrenzung als Schmalbandlautsprecher anzusehen ist, 
falls man von ihm kleine Verzerrungen bei grosser Leistung fordert. 


Über die Ausbreitung „exponential-begrenzter“ Wellen von endlicher Amplitude 


Als Einleitung sollen zunächst die bekannten Gleichungen der Luftbewegung in 
einem Horn angeführt werden, wobei angenommen wird, dass diese Gleichungen 
(für Wellen von infinitesimal kleinen Amplituden) durch Weglassen von als klein 
angesehenen nichtlinearen Gliedern bereits linearisiert wurden: 


У lun Du ci AS 
die Eulersche Gleichung: 00 3r fig eh 0, 
| дъ мез Alben ee 
йе e WE с ot tug ar о 
die thermodynamische “Gleichung: др =у до 
Po Qo 
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Von den drei auftretenden Unbekannten u, бр, до interessieren uns hier haupt- 
sächlich nur die zwei ersten, nämlich die Geschwindigkeit der Partikel und der akus- 


tische Druck. 
Nach erfolgter Trennung der Variabeln werden folgende Differentialgleichungen 
erhalten: 


92у (3 2) ди 9 (3 2 о 9u 


да? S Ox] 9x | 8x VS Ox] — yp, 08 
Эр С НЕЗ дА AE 
9x? S dx) ga урь ott | MA GS 


Fine bedeutende Erleichterung bei der Auflösung dieser Gleichungen bietet die 
Annahme, dass der Hornquerschnitt 5 exponential abnimmt, nämlich 


S = 50.6 
Dann ist 

ER 

S 8x q 


und die beiden Differentialgleichungen werden identisch 


dtu | du _ 18u 
0x33 49 сї 912 
923p 9др _ 1 Ә?др 


Эх? 4 Ox с ge 


Die partikulären Integrale dieser Gleichungen nehmen in reeller Schreibweise 


die Form an 
= 2U, e-«ax? = 2 — 2 
u = 1 
oe cos | o£ — x = 


2 


= = о g 
др = 2 P, е–!/2 cos КЕ EAR n +9). 
Sie entsprechen der Wellenausbreitung in Richtung der positiven Koordinatenachse 
x, die Phasenverschiebung ф hängt von den Randbedingungen ab. Falls wir die Kompo- 


nente berücksichtigen -wollen, die der Wellenfortpflanzung ins Innere des Hornes 
(d.h. in die Richtung abnehmender x) entspricht, dann haben wir 


— 
U = 2 Du e-#!2 cos (= = - z) dis 


2 Н 
+ 2 U,_ € 74*1? cos 7 Es 2 z) 
c 4 


nn 
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1. Einleitende Bemerkungen 


Bevor wir das eigentliche Problem anfassen, wollen wir noch zunächst kurz den 
Begriff der Exponentialbegrenzung erörtern. Eine Welle von endlicher Amplitude 
kann in Stromröhren zerlegt werden. Die Seitenbegrenzung des Flusses ist von mass- 
gebender Bedeutung für die Form dieser Röhren. Sie bildet nämlich eine Fläche, in 
welcher keine Quellen des Mediums vorhanden sind. Durch passende Wahl der 
Krümmung dieser Seitenbegrenzung kann man eine Welle von folgendem Aufbau 
erhalten: 

Zieht man Querflächen senkrecht zu den Stromröhren, so dass die Wellenphase 
in jedem Schnittpunkt jeder Röhre mit einer dieser Flächen dieselbe ist, und nimmt 
man dann die längs einer beliebigen Röhre gemessene Länge als Parameter dieser 
Phasenflächen an, dann soll die Begrenzung Flächen ausschneiden, deren Inhalt expo- 
nential mit dem Längenparameter der auserwählten Röhre wächst. 

Aus den noch im XIX Jahrhundert durchgeführten mathematischen Unter- 
suchungen von Hugoniot (1889) könnte man schliessen, dass die Exponentialbegren- 
zung eine der natürlichen Stromröhrengestalten bei grossen Störungen des Mediums 
ist. Die Gleichung für die Querschnittsfläche der Stromröhre bei idealer Flüssigkeit 
besitzt nämlich die Form: 

GAS А du u? 
Ge tS (1-4) =0 Gees 
Bei sehr grossen Amplituden kann 1 als sehr klein im Verhältnis zu u?/c? gestrichen 
werden. Dann erhalten wir nach erfolgter Integration der letzten Gleichung: 
E (2) 

Diese Begrenzung entspricht nicht genau dem Begriff, welcher in der Technik 
als Exponentialhorn bekannt ist. Im letzten Fall vergrössert sich nämlich der zur 
Achse senkrechte, flache Querschnitt exponential mit der Länge des Hornes; nur die 
sich langsam erweiternden Horne können als exponentiale Flussbegrenzung angesehen 
werden. s 

Das Problem der Wellenfortpflanzung in einer Exponentialbegrenzung ist.vor 
allem deswegen interessant, weil die Welle im schmalen Teil der Begrenzung die 
bewegliche Randbedingung erfüllen kann, was bei Wellen von endlicher Amplitude 
immer vorkommt. Gleichzeitig wird in dem breiten Teil die Welle unendlich klein, 
was dort die Anknüpfung an die Theorie der Infinitesimalwellen des Mediums erlaubt. 

Wie bekannt, bestehen in der Anwendung der nichtlinearen Theorie flacher 
Wellen grosse Schwierigkeiten bei der Betrachtung der Wechselwirkung der Welle 
mit irgendwelchen unbeweglichen Randbedingungsträgern, z. B. mit der Impedanz 
eines im Wege der Welle stehenden festen Hindernisses. Auf Grund der Nichtlinearität 
der Gleichungen ist es unmöglich (wie es in der linearen Theorie üblich ist) die ganze 
Welle aus einem hinstrómenden und einem zurückkehrenden Teil zusammenzusetzen. 
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Die Sache lässt sich dagegen leicht in der Exponentialbegrenzung durchführen, 
falls das oben erwähnte Hinderniss im breiten Teil, wo die Welle infinitesimal ist, 
eingesetzt wird. Die ausgestrahlte Welle unterliegt nämlich beim Fortschreiten im 
schmalen Teil der Begrenzung einer gewissen Entartung, ühnlich wie das bei den 
ebenen Wellen der Fall ist. Diese Entartung nimmt aber stündig ab, als die Welle nach 
dem breiten Begrenzungsteil fortschreitet. Die Welle wird dann in der festgesetzten, 
entarteten Form reflektiert und unterliegt beim Wiedereindringen ins Innere des 
Hornes einer erneuten weiteren Entartung, die im schmaleren Teil der Begrenzung 
grósser wird. 

Weiter muss folgendes berücksichtigt werden: um grosse Komplikationen und 
mathematische Schwierigkeiten beim Auflósen der Gleichungen auf Grund der Theorie 
der Wellen von endlicher Amplitude zu vermeiden, ist es ratsam, einen anderen Aus- 
gangspunkt als bei der linearen Theorie anzunehmen. Man gibt am besten die Glei- 
chungen in „substanziellen“ (Lagrange'schen) Koordinaten an, von denen man sagen 
kann, dass sie „mit dem Medium fliessen“. Wenn eine Partikel des Mediums (Kot- 
schin, Kibel und Rose 1948, SS. 48..., 352... ; Kutilin 1947, SS. 158...) die Koordinate 
x in der Ruhelage hatte, dann wird sie nach erfolgter Verschiebung mit derselben 
Zahl x — x bezeichnet. Die wahre Verschiebung £, in geometrischen Koordinaten 
gerechnet, ist eine Funktion der Koordinate x und der Zeit, sie beträgt also &(x,t). 

Die geometrische Lage der Partikel x ist im Augenblick t 


x = x + é(x,t) (3) 


2. Ausbreitung der Welle innerhalb einer Begrenzung von bekanntem Verlauf 


Die Wellenflüchen in der Begrenzung bilden eine Familie, die die Achse senkrecht 
durchschneidet (Abb. 1). Ein Teilchen des Mediums, das sich innerhalb der Begrenzung 
bewegt und in Abwesenheit einer Störung sich zwischen den Flächen S(x) und S(x + 
+ dx) befindet, verschiebt sich unter der Wirkung der akustischen Welle in eine neue 
Lage, die durch die Flächen S(x + &) = S(x) und S[x + & + d(x + El begrenzt 
ist. Die Dicke der verschobenen Schicht beträgt also dx = d(x + ё) = ( ЕРЕ 2j dx. 
Es sei nochmals ausdrücklich darauf hingewiesen, dass die Variable x die Koordinate 
der Anfangslage des Luftteilchens im Inneren der Begrenzung ist und dass sie von 
der Zeit unabhängig ist. Die Variable E dagegen ist eine Funktion von x und auch t 
und stellt die Verschiebung des Teilchens aus seiner ursprünglichen Lage unter dem 
Einfluss der akustischen Welle dar. Die Anfangsdichte sei durch о(х) = 00 = const 
bezeichnet, ähnlich sei р(х) = pg = const. Die Dichte des verschobenen Mediums 
о(х + EI ist mit dem Druck durch die adiabatische Formel verbunden (y = C,/C,) 


р(х) _ | sah 


Ser Joer (4) 
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Die Kontinuitätsgleichung kann > geschrieben werden 
S(x)e(x) = 5(х)(1 + Eo), (5) 
wo EI = 9&/9x die Deformation ist. 


Als dritte Gleichung wird diejenige des Gleichgewichts der Druck - und Trag- 
heitskräfte, die auf ein Element der Stromröhre wirken, angenommen: 


2% рб) 
aa 0 S 


Das System der Gleichungen (4), (5) und (6) mit den Unbekannten р, о und & beschreibt 
das akustische Feld von endlicher Amplitude innerhalb einer Begrenzung von bekann- 
tem Verlauf. Wir nehmen jetzt gewisse Änderungen vor, um die Gleichung der Ver- 
schiebungsfortpflanzung & zu finden. 


Abb. 1. Zur Geometrie der kleinen Verrückungen in Exponentialbegrenzung. 


Die Kontinuitätsgleichung kann gleich in einer solchen Form geschrieben werden, 
dass die bekannten Grössen oy und 5 auf einer Seite-zu stehen kommen: 


5(х) 

"ag ` 

Das letzte Glied der Bewegungsgleichung (6) lässt sich folgendermassen umformen 
(vgl. Earnshaw 1859) | 


= 0(х) 1+). 


Эр(х) _ Әр(х) Әх 


Ox уе суы. 


aber 
ae PZ IT Уве №1 
9х Ox ALE 148” 
Ox 9% 
also 
Р Әр(х) 9p), 1 (7) 
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Aus (6), (7) und der umgeschriebenen Kontinuitätsgleichung folgt 
92E S(x) Әр(х) 
OM — 0. 8 
92:3 Te Әх (8) 
Um die Ableitung zu berechnen bedient man sich der thermodynamischen Glei- 
chung (4). Löst man sie in Bezug auf p(x) und setzt aus der Kontinuitätsgleichung 


o(x)/o(x) ein, so bekommt man 


р(х) = zo 


S(x) Le 
= (1+ =’) 
ET 
Diese Grösse wird jetzt nach x differenziert, indem man berücksichtigt, dass & auch 
eine Funktion von x ist. Dann wird 


aso 2 
9p(x) _ = Уро 9х S(x) | 9x? 
дх [se Y e 1+ é 


Jetzt wird diese Grösse in (8) eingesetzt und уро wird vermittelst der bekannten Be- 
ziehung уро/оо = c?, wo c die Schallgeschwindigkeit für Wellen mit unendlich kleiner 
Amplitude ist, durch бес? ersetzt. So findet man die Ausbreitungsgleichung für die 
Verschiebung & in einer beliebigen Begrenzung: 


2 2 
1 92 1 9х S(x) 9x? 2 
S(x) S(x) 


3. Die Exponentialbegrenzung 


Gleichung (9), die auf den: ersten Blick ziemlich unübersichtlich erscheint, ver- 
einfacht sich bedeutend für die Exponentialbegrenzung, für welche gilt: 5 = 50е“. 
Laut der ursprünglichen Bezeichnungsweise setzen wir jetzt wieder x + & an Stelle 
von х ein. Die Grössen S(x + §)/S(x) und (9/9х)(5( + £)/S(x)) werden jetzt durch 


die Entwicklungen 


Se Ө ат urn 

S(x) 
9 SEHE) _ Ê y op 
SRS agnor | reae +. 

ausgedriickt und die Ausbreitungsgleichung nimmt folgende Form an 

> de 

1 2x? 2 

Е x 192€ PS (10) 


| 222 SH Ela s BEE) 
D + qê + uw + >) ИД wu EB 
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4. Integration der Gleichung der Verschiebungsfortpflanzung 


Diese Integration wird nach der Methode der sukzessiven Annäherungen aus- 
geführt, indem wir von der Randbedingung für x = 0 ausgehen. Im Punkte x = 0 
erzeugt die Kolbenmembran eine akustische Welle, indem sie dort die folgende 
Verschiebung in der Luft bewirkt 


оа еіп == К A sin wt (11) 


A = beliebige Amplitude in dimensionsloser Schreibweise, wobei das Mehrfache der 
Wellenlänge 2/2 = {Го die Längeneinheit für die dimensionslose Amplitude bildet. 
Wenn die sukzessiven Annäherungen der Funktion &(x, t) eine Potenzreihe mit 


der Variable A: 


Et) = — [Афу + Ap, + Абру +...) wor ф = фы) . (12) 
bilden sollen, dann müssen für x = 0 die Funktionen ф folgende Anfangsbedingungen 
erfüllen 

ф1(0,1) = sin wt 
(0, = 0 
2(0,t) = 0 


++ ж ж о э» э о е 


Es wird auch angenommen, dass am anderen Ende keine Reflexion stattfindet, so dass 
die Funktionen sich in positiver Richtung der x— Achse fortpflanzen. 
Das Einsetzen der zwei ersten Glieder der Reihe (12) in die Gleichung (10) ergibt 


1],9?p | „Фа ЗА? U, КЕС m 
L4 3n Lu anf er orte obe 

OUT 99i 2992 cA Эф, cA dpa _ с2А? 9g, Эф; | 

а (4B +4 DE |e а | айл» 
Diese Gleichung lässt sich in zwei simultane Gleichungen spalten, von denen die erste 
A, die zweite А? enthält. Jede der beiden Gleichungen F(x,t) = 0 und G(x,t) = 0 
muss einzeln erfüllt sein, nur dann kann AF(x,t) + A?G(x,t) für beliebige A ver- 
schwinden. Die erste Gleichung enthält nur die Funktion e: 


| ox? gx OR gum e 
Sie ist uns schon aus der Exponentialhorntheorie bekannt. Das Integral dieser Glei- 
chung, das die Anfangsbedingung erfüllt, lautet 


5E pi 
sin (o —% e T SL (14) 


p m _ 19%, — аз) 


і و‎ 
p(t) =e ? 4 
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Die zweite Gleichung, deren Glieder A? enthalten, ist 


9%», дф _ 199 сдф 9% | q а, 1 رو‎ d az 
3m LEE 5852/0: 0: a кайыл ома 


Sie ähnelt einer Gleichung von bekanntem Typus, enthält aber ein inhomogenes 


Glied (das mit y bezeichnet werden soll), das eine periodische Funktion von der 
Frequenz 2w ist. Man ersieht daraus, dass dieses inhomogene Glied aus Produkten 
zweiten Grades bekannter periodischer Funktionen der Frequenz œw aufgebaut ist. 
Wir rechnen das Glied aus und erhalten nach entsprechenden Kürzungen 


d ay —1)® ر‎ gy = 5) o* + ge? “e082 (ot - 3 o? a 


Vili Si 4c Acw "ET 


c? 4 


(y + Џо? — q?c? w? q? ‹ = б g 
НЕ I 2 AR sin 2 | wt — x =) 4 


Das Integral der Gleichung (15) besteht aus zwei Bestandteilen. Der erste Teil ist 
das Integral einer homogenen Gleichung (Funktion vom Typus (14)) von beliebiger 
Frequenz. Um die Anfangsbedingung zu erfüllen, wird ihm die Frequenz Zo zuge- 
schrieben. In Analogie mit den früheren Gleichungsintegralen wird er folgendermassen 
geschrieben: 


eh E ER 
gue i | cs (2o Sg s -5) -- D cos (2o — «V 5 -$Jjae 


Die Gestalt des zweiten Bestandteils wird auf die des inhomogenen Gliedes zurückge- 
führt: 


x 2 g? pel 
фаз = vele: Q sin 2 [o — Ex е, + Rcos2 (or ТЕ) 
A 


(17) 


Nach dem Einsetzen von фу» in die Gleichung (15) haben wir: 


№: A 
Ws Se ie AL a "ET: дус _ (у+0о 
Paa EU E (Epp dps "ri age = (18) 


c 
Man sieht daraus, dass für x = 0 das Glied 554-0, seine Koeffizienten jedoch genau 
bestimmt sind. Es folgt daraus, dass das Erfüllen der Anfangsbedingung ¢,(0,t) = 0 


durch die Summe 9, = Pa + Pag nur für entsprechend erwählte Koeffizientwerte 
im Integral Ф, möglich ist. Es folgt, dass 


C=-0; D=-R (19) 


5. Intensität der zweiten Harmonischen in der Exponentialbegrenzung 


Es ist interessant zu erforschen, wie gross die Intensität der zweiten Harmonischen 
am Ende der Exponentialbegrenzung ist. Sie wird durch das Verhältnis der Quadrate 
der Geschwindigkeitsamplituden des Teilchens beim Ausgang gegeben. Mit Rück- 
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sicht darauf, dass die Frequenzen p, und фі entsprechend 20 und о betragen, ist 


dieses Verhältnis 
le Aaron 44° фо]? 
Пё rd Uri is oA э Ж m 
фл op, | xa KAR /z=ı 


wo |ф| die Amplitude der sinusoidalen Grösse ф bedeutet. Wir schreiben in expliziter 


Form 

_ a 2 2 ql 
Oe Sar sin (or — | | / 3-2) CET 
| polz=1 = | Par + Фа = 


| 


-g 
= ] / 0 / 4а? д? 
do (2-1 D -E)- Ros (201 - 1 I al | 
AD 2 
ev [PL + Qin 2 o - (ee -£j 4 Roos? (= VEH 
@ 4 ce 4 


Nach einfachen Umrechnungen erhalten wir 


41 gl 


lpslz-;— |e “PI| +e ?B—e 2-B, 


wo 


" 
B? = (0° + R3) d Legd — 2e 2 cos (x — al (21) 


У 1 ЕЕ en =. (22) 


Die relative Intensität der zweiten Harmonischen beträgt also 
n = (24B)? - (23) 


Man sieht daraus, dass diese Intensität mit dem Quadrate der dimensionslosen 
Amplitude der erregenden Kolbenmembran wächst. Für eine gegebene Amplitude 
ist sie der Grösse B? proportional, welche eine Funktion der Länge / der Begrenzung, 
des Exponenten q und der Frequenz w ist. 

Wir führen die Bezeichnung Ё = w/w, ein und berücksichtigen, dass die kritische 
Frequenz wg = qc[2 ist. Dann haben wir nach einfachen Umformungen 


gl 
B = [оо k? — E + Sek £ Leg де SC (E al (24) 


und 


16 
1 
Für I> A) = 4r/q sind die Glieder e? und e~ viel kleiner als 1; man kann sie 


also völlig vernachlässigen und В nur als Funktion von k behandeln. Abb. 2 stellt 


444 EH RY Ве stat 


В? als Funktion von k dar. Diese Funktion ist stets positiv und erreicht ihr Minimum 
für А = (10/9)"*. ) 

Selbstverständlich um vom Koeffizienten В? zur wirklichen Leistung überzu- 
gehen, muss man im vorliegenden Fall die aus dem Diagramm abgelesene Grösse 
mit 4А? multiplizieren, wo A die dimensionslose Ablenkungsamplitude ist: 


E 


c 


und |&| die Ablenkungsamplitude der Kolbenmembran. 

In der Exponentialbegrenzung breitet sich die Welle mit grosser Gestaltsentar- 
tung aus in dem Bereich der Frequenzen, die niedriger als die Abschneidefrequenz 
sind. Für die Abschneidefrequenz beträgt der Koeffizient ungefähr 0,003, im Inter- 
valle 1<k<2 erreicht er sein Minimum. 


В? 


Abb. 2. Der dimensionslose Faktor der Amplitude der zweiten Harmonischen in Exponentialbegrenzung 
als Funktion der dimensionslosen Frequenz. 


Um die Exponentialbegrenzung zum Bau eines Grossleistungslautsprechers an- 

zuwenden, dürfen als Erregungsfrequenzen nur diejenigen Frequenzen zugelassen 
werden, die sich im Intervalle befinden, für welches B? klein ist. 
Insbesondere, obwohl der Wirkungsgrad des Lautsprechers unter der Abschneide- 
frequenz sehr klein ist (d. h. die Ausstrahlungsbedingungen sehr ungünstig sind), 
werden die Verzerrungen so gross, dass man die Ströme niederer Frequenz nicht an 
die Erregermembran zulassen darf. 

Es zeigt sich also, dass ein Lautsprecher, der eine nach dem Prinzip der Exponen- 
tialbegrenzung hergestellte Welle ausstrahlen soll, ein Schmalbandlautsprecher ist, 
wenn man ihm die Bedingung kleiner nichtlinearer Verzerurngen bei grossen Leistun- 
gen auferlegt. 

Dem Leiter des Instituts Herrn Prof. Dr S. Szczeniowski danke ich herzlichst für 
die Schaffung der. Bedingungen, welche die Entstehung dieser Arbeit ermöglicht 


haben, ebenso Herrn Prof. Dr M. Kwiek für dauernde Unterstützung und kritische 
Diskussionen. 


le ы 
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КРАТКОЕ СОДЕРЖАНИЕ 


Г. Рифферт, РАСПРОСТРАНЕНИЕ ВОЛН С ОКОНЧАТЕЛЬНОЙ АМПЛИ- 
ТУДОЙ В ЭКСПОНЕНЦИАЛЬНОМ ОГРАНИЧЕНИИ. 


В настоящей работе было изунено распространение акустической волны 
в экспоненциальном ограничении. С помощью криволинейных координат этот 
вопрос удается сформулировать, как одноизмерную задану. Волны образуются 
под действием гармоничных колебаний поршня, находящегося в узкой насти 
ограниченности. 

Исходную тонку составляют представленные в координатах гидродинамине- 
ские управления Лягранжа, решение которых основывается на развитии всех 


2n|5 
величин согласно степени относительной амплитуды A = ai ‚ где él является 


амплитудой гармонинеских колебаний поршня, 1 — длиной волны, которая со- 
ответствует круговой настоте w колебаний поршня. 
В первом и втором приближениях — образующиеся гармонические состав- 
ные волны — имеют частоту w или 20). 
59% 
Волна с настотой равной о угасает с функцией c 2 вместо того волны 


NE 
C настотой 2% частично угасают C функциейе 2, а настинно C eg, 


В публикации указаны формулы для выснитания соотношения интенсивно- 
сти обеих волн при унёте ограничения. Это соотношение оказалось пропорцио- 
нально к квадрату относительной амплитуды. 

Окончательным выводом является то, что громкоговорителем с экспонен- 
циальным ограничением, от которого требуется небольших искажений при боль- 
шой мощности, нужно снитать громкоговоритель узкополосовый. 
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THE YIELD OF ANTI-STOKES FLUORESCENCE 
OF VERY VISCOUS DYE SOLUTIONS 


By REGINA DRABENT AND Danuta FRACKOWIAK 
Physics Department, Nicholas Copernicus University, Torun 
(Received January 23, 1955) 


The results of measurement of the relative yield of fluorescence of fluorescein in 
glycerol and glucose are given. In the first case the yield decreases with increasing wave- 
length of the exciting light, in the second case, however, it does not — remaining practi- 
cally constant. An explanation of the experimental results is put forward, according to 
which the course of the yield as a function of the exciting wave-length depends on the 
relative concentration of non-luminescent dimers. 


Introduction 


The dependence of the fluorescence efficiency on the wave-length of the exciting 
light was investigated in detail by S. I. Vavilov. According to this author the energy 
yield of fluorescence (i. e. the ratio of the emitted to the absorbed energy) of aqueous 
solutions of fluorescein increases linearly with the wavelength of the exciting light 
up to the anti-Stokes region, where the efficiency drops abruptly. This steep drop 
of efficiency was confirmed for different solutions by several investigators. Vavilov 
has expressed the opinion that quite generally: р 

(1) the energy yield of luminescence can never exceed the value of 1 and 

(2) the lower the temperature of the solution the steeper must be the drop of 
efficiency in the anti-Stokes region. 

Vavilov (1945, 1946) and Adirovich (1945) have argued that the above statements 
result from the second law of thermodynamic if the quantal character of the pheno- 
menon is taken into account. Several attempts were made to explain the phenomenon 
from molecular point of view. Pringsheim (1949, p. 308) has put forward an expla- 
nation, according to which the drop is due to the perturbation by the interaction of 
the excited molecules with the adjacent molecules, which is particularly strong in cases 
of absorption in wings of the absorption band, and thus in the case of absorption in 
the anti-Stokes region. Vavilov (1927) and Jablonski (1931) suggested that the drop 
of the yield at greater wave lengths of the exciting light may be due to a very weak 
“inactive“ absorption band. Later Jabłoński (1933, 1935) expressed the supposition 


(447) 
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that this band may be possibly due to a forbidden transition from the ground state N 
to the metastable state M. The last explanation was criticised by Pringsheim (1949), 
who pointed out that the “transition from N to M has a probability almost 108 times 
smaller than that of transition from the ground state to the fluorescent state F.“ The 
M—N band would thus be much too weak to account for the observed drop. Although, 
according to Jablonski, in the case of solutions of low viscosity the probability of N—M 
transitions may be enhanced very considerably by rather strong perturbations of lu- 
minescent molecules by adjacent molecules of the solvent, the perturbations being 
strong enough to quench phosphorescence, he has lately dropped his former explana- 
tion and put forward a new one. Most of the luminescent molecules polymerize. In many 
cases the polymers, say dimers, are not luminescent. The absorption of light in an 
absorption band belonging to such dimers does not lead to luminescence, it is „inac- 
tive^. The maxima of absorption bands of dimers, in the cases quoted by Pringsheim 
(1949, p. 354) and Levshin (1951, p. 177) are situated on the short waves side of the 
maximum of the absorption band of monomers (“active“ absorption band), but the 
dimer absorption band is much broader than that of monomers — its long wave tail 
is superimposed on that of monomers and even extends beyond it. Thus the absorption 
of light in the region of the long wave tail (the anti-Stokes region) may be very ineffi- 
cient indeed. The relative concentration of dimers depends on concentration of the 
solute and thus the shape of the “absorption band“, being actually a superposition 
of bands of monomers and dimers, depends on the concentration of the solute, while 
that of the fluorescence bands does not depend on it considerably. At higher tempera- 
tures the dimers dissociate into monomers. Hence the decrease of yield becomes less 
steep at higher temperatures. These are many papers concerned with experimental 
investigations of the relative yield of anti-Stokes fluorescence of dye solutions.! 

Apart from papers of Vavilov (1927) and Jablonski (1933), who was the first to 
prove Vavilov's results by means of an objective method, a great number of papers 
were published in which similar results were reported for many different dye solutions 
(cf. e. g. Solomin 1941). 

Originally the aim of the present work was to find out whether the N—M transition 
can be made responsible for the decrease of yield in the anti-Stokes region (the experim- 
ental results concerning the influence of temperature on the slope of the yield curve 
were unknown to the writer at the time when this work was started.) It was supposed 
that, if the inactive band is due to N—M transition, it may be weaker in very viscous 
phosphorescent solutions than in non phosphorescent solutions, in which stronger 
perturbations may enhance the transition probability. First the fluorescein solution 
in boric acid was tried. Such solutions show a very intense phosphorescence, but soon 
become turbid and thus unsuitable for investigations, and had to be abandoned. We 
then investigated a solution of uranin in anhydrous glycerol. This solution, however, 
did not show any phosphorescence, at all. Since we suspected that the lack of phospho- 
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1 For earlier papers cf. Pringsheim (1949, p. 208). 
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rescence may be due to quenching of excited molecules by oxygen, we tried to disozi- 
date the solution, but did not obtain any change in properties of the latter. The relative 
energy yield of fluorescence 7/7 as a function of the frequency > of exciting light 
appeared to decrease according to Vavilov’s expression 

N = en 

Л = eb”) | (1) 
However, the slope of this function appeared to be less steep, than that to be expected 
on ground of the assumption 6 =h/kT, this behaviour being similar to that observed in 
other solutions so far investigated. 

We next investigated fluorescein solution in sugar (saccharose) as an example 
of a very viscous, glassy solution. This solution did not phosphoresce, probably because 
of the rather poor quality of the used fluorescein (commercial fluorescein, provenance 
unknown). There was no drop of the yield observed for exciting wave-lengths between 
5000 and 5450 А. U. The yield did even increase in the anti-Stokes region, although 
slower than in that of the Stokes excitation. We tried to go further into the anti-Stokes 
region with exciting wave-lengths (according to Pringsheim the peak of fluorescence 
band lies at 5270 А. U. thus excitations with any 2 > 5270 А. U. constitute the anti- 
-Stokes excitation) hoping to find the yield drop in a farther region. But the necessity 
of using in the experiments long exposure times, during which the glassy solution 
began to recrystallize, prevented us from obtaining reliable results. The primary light 
scattered in the turbid solution falsified the experimental results, which thus had to 
be rejected. 

Further experiments were made with fluorescein solutions in glassy glucose. 
These solutions proved to be phosphorescent and rather stable. It appeared that by 
excitation by different wave-lengths varying from 5200 up to 5760 À. U. the energy 
yield of fluorescence remained practically constant (see Fig. 5). This result is to be 
expected if the drop of the yield is due to the influence of the absorption of light by 
non-luminescent dimers. ў 

At higher temperatures, at which the glucose solution is prepared, the relative 
concentration of dimers becomes very low. The cooling proceeds rather quickly. At 
lower temperatures the enormous viscosity preverits the solution from reaching a new 
equilibrium state — the relative concentration of dimers corresponds to the equili- 
brium at the temperature to which the solution was heated (about 160°C) before the 


cooling began.” 


Experimental 


A photographic method, first applied by Jabłoński (1933), was used. Fig. 1 shows 
the arrangement. The image of the light source (a low voltage incandescent lamp fed 
by a battery of accumulators with a better voltage constancy than + 0,01 V) was 
EE RR E ав 


2 The same applies, of course, to dye solutions in saccharose, 
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projected on the slit of a double Hilger-Muller monochromator. The light beam 
leaving the exit slit of the monochromator M is made parallel by a set of additional 
slits S,, 5», and then enters the fluorescent solution. It there produces a decaying 
(because of the gradual attenuation of the primary light by absorption) fluorescent 
beam, which was photographed by means of a ,,Contax“ camera with its axis perpen- 
dicular to the direction of the primary light. In order to get larger photographs than 
it is usually possible with the above camera, 
the tube ofthe camera has been lengthen- 
ed by insertion of an additional piece of 
tube. “Ultra-Rapid“ Agfa films have been 
used. 

The intensity of the beams of 


light of various wave-lengths given by 
A the monochromator was measured by 
means of a Moll thermopile with a 
Moll galvanometer for different widths 
of the monochromator slits, and extra- 
polated to the width actually used in 
C experiments. 
A In orderto obtain the sensibility curves 


of the films every beam was photographed 
Fig. 1. Experimental arrangement. B — bulb, > ү Se 


І — lens, $,, S, — slits, M — monochromator, Е : 
С — container with the fluorescent solution, reduced in one of the exposures in a known 


А — photographic camera,  — wire net. proportion by interposing a wire net 
before the lens L. From such a pair of 
photographs the sensibility curve of the film can be obtained (cf. Jablonski 1933). 
From a set of photographs made on a single film one gets all the data needed for calcu- _ 
lation of the relative yield of fluorescence 7/7, provided the relative intensity of the 
primary light /(А;)/1(40) for various wave-lengths A is known (this relative intensity, 
as mentioned above, was measured by means of a thermopile, the subscript “о“ denotes ` 
the quantities belonging to the comparison yield nọ, the subscript *1* — those concern- 
ing the yield to be compared). These data are: the intensity of the fluorescent beams 
F(A,x) as function of the distance x from the entrance of the primary light I(/) into 
the fluorescent solution for different exciting wave-lengths A, and the corresponding 
absorption coefficients u(A). The relative yield is to be calculated by means of the fol- 
lowing expression (cf. Jablonski 1933) f 


11 Es Е (21 EA I, (Ao) Ho (Ao) еш): Ho(Ao):x0 (2) 


no Fo (Aosxo) А (А) (Ar) 


Results 


twice, the exciting light intensity being 


(a) Uranin in glycerol 
The glycerol was dehydrated so as to contain 2 per cent water only. The concen- 


M E» 


The Yield of Anti-Stokes Flourescence 451 


tration of uranin was 1.6Х 10-4 g cm^?. The widths of the monochromator slits were 
0.25 mm. The corresponding width of the spectral region transmitted by the mono- 
chromator was in our conditions about 50 A. U. Fig. 2 shows the absorption coefficient 
of the investigated solution versus wave-length. 

The relative energy yield 7/7, (the ratio of 
yield with excitation by different wave-lengths 
to that by excitation by 7) = 5100 A. U.) is 
given in Table I and shown in Fig. 3, where 
also curves representing logarithms of the yields 
as function of /, obtained by other investigators 
for aqueous solutions of fluorescein are given 
for comparison. The experimental curve along 
with that calculated by means of Vavilov’s ex- 
pression (1) is given for comparison in Fig. 4 5000 500 500 5600 S800 
(expression (1) requires Log (n/n,) to be linear 
function of v— »,). 


Fig. 2. u versus A 
© for uranin in glycerol (1.6 x 10-4 = cm 3) 
x for fluorescein in glucose (1.45 x 10-4 g/g). 
A for fluorescein in glucose (1.43 x 10-3 g/g). 


(b) Fluorescein in glucose 

Fluorescein sodium (of Coleman Bell Co 
Norwood О, U.S.A.), soluble in water, was used. The method of producing glucose 
solutions was that of Starkiewicz (1929) for sugar solutions. The concentrations used 
were: 1.45 x 10-4 g of fluorescein 
per gram of glucose (table II set a) 
and 1.43 x 10-4 g/g (table II set b). 
The spectral width of the exciting 
light was, as before, about 50 А. U. 
The fluorescence did not contain 
any perceptible trace of scattered 
primary light, as was checked by 
special experiments. The position 
of the fluorescence band as well 
as its shape do not change if the 
fluorescence is excited by light fil- 
tered by interference filters with 
maxima of transmission at 5172 


Ä.U. and 5875 A. U. respectively. 


-10 
9000 5200 2400 


Fig. 3. Log of relative yield of fluorescence = Log (77/70) Conclusions 
as function of wave length A of exciting light. 


The drop of yield of fluores- 
cence usually observed in the anti- 
Stokes region could not be ob- 


© fluorescein in dehydrated glycerol (present experiments)’ 
x aqueous solutions of fluorescein (Jabłoński), 

7 aqueous solutions of fluorescein (Vavilov), ; 
А aqueous solutions of fluorescein (Szezeniowski). served in the case of fluorescein 
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in glucose. This result can be explained by assuming that the drop in question is due 
to the preponderance of absorption by non-luminescent dimers over that of lumines- 
cent monomers in the anti-Stokes region. The glassy solution is prepared by boiling 
the starting product at a temperature of 160°C, at which most of the dimers “disso- 
ciate“ into monomers. The concentration of dimers remains practically unaltered 
when the solution is cooled, because of its great viscosity. It seems desirable to 


5400 5600 5800 
Fig. 4 Fig. 5 
Logn/no versus frequency of exciting light. The relative yield 7/7, as a function of the exciting wave- 
© fluorescein in glycerol, length A for fluorescein in glucose: X set а, A set b 


A X fluorescein in glucose, 
Dashed curve — theoretical curve 7/70 = 
= ej», where В = НКТ. 


investigate the influence of the boiling temperature used during preparation of the 
solution on the yield curve. 


The explanation of the drop of yield in the anti-Stokes region as due to the absorp- 
tion by non-luminescent dimers (Jablonski 1954) seems to be more plausible than 
that by forbidden N—M transitions (Jabłoński 1933 — The last explanation was 
criticized by Pringsheim 1949 and lately rejected also by Jabłoński 1954). 


The present results are not very precise (the experimental errors are different 
for different wave-lengths of the primary light and may be as high as 20 to 30 per 
cent), nevertheless the general trend of the dependence of the yield on the exciting 
wave-lengths remains in all sets of experiments (which were repeated several times) 
the same. Further experiments (are being carried out. 


We wish to express our sincere gratitude to Professor A. Jablonski for suggesting 
the problem and for giving experimental and theoretical help. 
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P. Драбент и Д. Фронцковяк, ВЫХОД ФЛУОРЕСЦЕНЦИИ ОНЕНЬ 
ВЯЗКИХ РАСТВОРОВ КРАСИТЕЛЕЙ В АНТИСТОКСОВОЙ ОБЛАСТИ. 


Измерен относительный выход флуоресценции раствора уранина в глицерине 
и флуоресцеина в глюкозе в антистоксовой области. 

В первом случае энергетический выход падает с возрастанием длины волны 
возбуждающего света, во втором в исследуемой области выход практически по- 
CTOAHEH. 

Эти результаты истолкованы Ha почве выдвинутого Яблонским (1954) объас- 
нения падения выхода/ Согласно с этим объяснением падание выхода вызвано 
наложением абсорпционной полосы нелюминенцирующих димеров. 
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In § 1 some general remarks are given on the Cauchy problem for the system: scalar 
meson field + particles. In § 2 an exact solution of the two-body problem in the case of 
a circular motion is given foru = 0. This solution is compared with the solutions in the 
cases of some different approximations. In § 3 some remarks are given on the rectilinear 
motion concerning the problem of the repulsive forces for small distances between the 
nucleons in the scalar field theory. It is shown that some physical contradictions result 
from the theory neglecting retardation. 


l. Introduction and General Discussion 


We shall start with the variational principle for n-bodies (nucleons) with the 
world lines &,(s,), s; being the proper time of the “i“-th body interacting by the scalar 
meson field p(x) corresponding to the meson of the mass m, (we shall denote u = m,c/h.) 

The action functional dependent on the meson field and the world lines of the 


particles of masses т; and charges g; enclosed within two space-like surfaces c, 


т 


and ос, Wo,0, is 


o 


Oa 2 i E 
W = Y [- «f ds; + gile f p (£) is) M | dy х Zap? — PP») (L1) 


а 


i 
The variational principle „ôW = 0,дф and û being equal to zero on o, and oy” 
yields the equations of field and motion: 


+00 

(7—2) (0 = An » gi if ds; ô, (x — 8), (1.2) 

d E — gile OY D = &lc Фа (8); iln (1.5) 
(455) 
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From Eq. (1.2) — (1.3) the conservation principles may be easily deduced. Let 
the energy-impulse tensor be 
1 1 
Ta = Е Фе — 5 (Ps, — we (1.4) 


4лс 


Ge 
D 


If T? vanishes at infinity & being everywhere finite, we have for the expression for the 
mpulse of the system „particles-field“ the invariant form 


Р, [o] = ee E| + f dese (1.5) 


i 
and for arbitrary о and x the conservation law is valid in the form 


óP.[o] _ 
So) — 


The system of Eqs. (1.2)—(1.3) is a „hyperbolic“ one, since the initial conditions 


(1.6) 


for the field and the world lines of the particles on an arbitrary су, e. g. ф(х) and 
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future motions of the particles. In fact, let A(x) be the function well known from the 

quantum field theory defined by: xj = 0>A(&) == 0, A(x)» = 6,(x); (O — и?) 

ХА O 


Let G(3)[o,] be the Green function parametrically dependent on o,, defined by! 


(C142) С (3)[oy] = — 4пб(х) 


The solution of the equation (07—02) U(x) = — 4ло(х) may be then expressed in 
terms of the initial conditions for U(x) on o, as follows: 


Ue) = | د-6 ریه‎ lode) + | dert) |2862) шуу 200) A زر‎ 
9y 
& | Й (1.7) 


1 
Thus, on solving Eq. (1.2) for ф(х) we obtain ф(х) expressed in terms of £,(s;) and the 
оси conditions for the field. On substituting this gx) into (1.3), we obtain for 


Es) a rather complicated system of integrodifferential equations with second order 
derivatives and kernels expressed in terms of the initial conditions for ф(х). 


1 For „flat“ og: % = op G defined d —p?) G (x — y) = — 470, (x — y); x? ао бы 
BC ў =~ eyr Ey 
= 0, AEP is given explicitly by: G = GG edic = A (x — y). 
"os 


The generalization for arbitrary o, is straightforward. 
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For the divergent coefficients of the 2nd derivatives usual mass renormalization 
procedure may be applied. The renormalized system of equations may be then reduced 
with the help of simple integrations to the system of pure integral equations containing 
parametrically the initial conditions for the world lines. 

It would be interesting to give a mathematical proof of the uniqueness and 


continuity of the solutions of these equations considered as functions of the initial 
i 


conditions for &, and ф(х). Probably the proof might be constructed with the help 
of the fix-point theorem in the conventional way. However, physically, it seems evident 
that the mentioned uniqueness exists.? From the physical point of view the „Cauchy 
problem“, (corresponding rather to scattering problems) .is not the suitable treat- 
ment for the case of interacting particles undergoing stationary motions. For the sta- 
tement of this problem the initial conditions for p(x) should be given, which for a sta- 
tionary problem are, however, determined by the motion of the particles in the whole 
interval of time. 

The problem is sometimes attacked on the basis of the radiation condition by 
the choice of the retarded potentials as the solutions of Eq. (1.2). Thus, evidently, the 
initial conditions for the field are simultaneously determined. The choice of the re- 
tarded potentials determines the radiation damping terms disturbing the stationary 
character of the motion. 

It is well known from the mathematical and physical point of view that this 
difficulty disappears in the case of symmetrical potentials (solutions of Eq. (1.2)). 

The aim of the present paper is, first of all, to investigate the influence of the 
peculiar transformation properties of the scalar field and the finite propagation time 


1 ə? | 

of the interaction (the term REC in Eq. (1.2)) on the character of the stationary 
с 

two body motions, as compared with the respective classical motions (с = оо). 


The effects of the peculiar transformation properties of the scalar field ф(х) 
(its appearance in the inertial term in the equations of motion), will be called the “kine- 
matical relativity“ (briefly “К. r.“). The effects connected with the finite propagation 
time of interaction will be called the “field relativity“ (briefly „f. r.“). 


+ 


2. The Circular Motions 


The rigorous solution of the problem of two-bodies of equal mechanical masses 
is possible in the case of circular motion for u = 0. At first we shall give the exact 
solution in this case and then the approximate solutions with neglect of the different 
terms connected with “К. r.“ and “f. r.“. These results are compared with the exact 


solution. 
REN dM Lc qo oe ЧЕ а E Ni ee a en 


2 J. Rzewuski (1953) proved the equivalence of a certain class of integrodifferential equations of 
second order with pure differential equations of the same order, proving simultaneously the uniqueness 
of the solutions of the integrodifferential equations by the same initial conditions. The problem discussed 
above is very similar to the problem solved by J. Rzewuski, it seems, however, to be mathematically more 


difficult. 
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(a) The case of retaining both “Е. т.“ and “У. г.“ foru = 0 (the rigorous solu- 
tion). The calculations for u = 0 seem to be of a rather academical character; they 
have, however, a certain meaning for the following reasons: 

(2) the calculations may be completely carried out and lead to compact formulae, 

(1) for small r we have e” —1 in the static Yukawa potential and it is hoped that 
the effects connected with и=Е0 are less important. 

Our system of the field equations and the equations of motion is 

+оо +оо 


Оф(х) = — Ang fe ô,(x — Ё) = дав [| 15, да(х — 5 (2.1a) 
d i = de 1 
as (тс — 51с$(5)) = = — g|cp,a($) (2.1b) 
2 
d 2 j 2 
ds; (me — ale) Ge = — aleve) (2.16) 


Since our problem is stationary, we choose for Green's function *G(x) = d(x, x’) 
After the mass renormalization procedure, (p(x) and ф, (x) should be substituted 
into (2.1b-c) respectively in the form 


+00 
19 (x) mM Гао == Zu (% Se EJ]; 19«(x) 
Zeta 
=g | EE — &) 2, éi (2.2) 
er | ee нй 
E 1 1 1 
=. | ds, 9 fpe ~ ©) (s — &)] +2 (ta Ё) 


—co 


Let us assume the motion to be circular, its velocity being v = ße: 
1 2 2 
ЧЕТ. COS e ERAS nee T2 Соз BEE En ^ 
үї = Gr Vl — в г 


1 
kee 


1 2 

©з = 0 > E,=0 , 
1 2 

Ee 51 GES 52 
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Wersee that 9, , = (È — &) (Е = £) — ү? ( SUSI Ай Cos B(s, — з») 
"ET И 
(2 


To verify the possibility of the postulated motion, we calculate 


Too +00 


Es) S A dsa д [012] = g V1 — | dvó[r*(v? — cos? fy)] 


+00 


= ҮГ BF йт | аре — оны 
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The second integral vanishes, since its integrand is the product of an odd and 


an even function e SUE: 
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On substituting the respective quantities үф and ф,, into (2.1 b), we find, 
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The first integral may be transformed into the final form 


where: 
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It is easily seen that the last equation is satisfield if 
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Similarly one can show that the relation (2.6 a) is sufficient to satisfy all the equa- 
tions (2.1 b,c) for our circular motion (2.3). 

For the sake of further discussion we shall introduce the dimensionless quantities 
g тС g? am, 52 + 274 
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= 0.01369. In this notation (2.6 а) may be rewritten as 
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To find f? = f*(x) we must calculate ДВ?). We have: 
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On introducing u = v—cos v, we find: 
dv 1 1 
dy эң ОЕ pone n 7 
Der, Dec l + Bsin vo 7o S 


where уу satisfies? vo — pcos vg = 0. 
However, we can obtain graphically В? = ß?(x) already from (2.6a) and (2.7). 
On taking уу as parameter we can write our curve in the following parametrical form 
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The curve (2.8) gives the velocity as function of the radius for our problem. 
The energy of the system: „particles + field“ in the general case may be shown to be 
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The last equation is valid for (2.1) in the case of the symmetrical potentials after re- 
normalization. It is obtained after substitution of the field expressed in terms of the 
world lines into the energy-impulse tensor and after performing the space integrations. 
With the help of the above equations the energy in the case of our circular motion 
may be calculated). 

It is interesting 10 compare this curve with the respective curves corresponding 
to some different theoretical approximations. We shall discuss the following cases: 
(b) The non-relativistic case (с = со in the field equations and the equations of 
motion). We obtain for u = 0 in the equation of field (the Coulomb case): 


В? = 20 | (2.9) 
and for p #0: 
В? = E e* (1 + x) (2.10) 


(c) The case of the potential being the zero component of a vector with retaining 
of »k. т.“ and neglect of the interaction of the magnetic type*. For our circular motion 


we obtain: 
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the expressions for the interaction forces. 
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4 We neglect Е LE “ar in the field equations and the contributions from the vector potential in 
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(d) The “k.r.“ effects in the case of the additional neglect of relativistic terms 
in the sources of the scalar field. It seems at first that the discussion of the effects 
of “k./r.“ will be reasonable, the limit coo being “consequently“ taken into account 
in the field equations. Eq. (2.1а) may be rewritten, as follows 
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On applying “consequently“ the limit coo in this equation, we find 
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On using for p(x) in the equations of motion the solution of the last equation and 
on neglecting the selfaction (the renormalization procedure may be applied), we find 
for our circular motion 
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Since — g? == > 0, the last equation may be satisfied only for r such that 


ra 
2 дг 
тс? — ef: > 0. For u = 0 we find: 

% 


Eq. (2.14) is reasonable for 2xyX:x-— oo. In the region OŞx<2x, our circular 
motion would be impossible. 


For u-F0 we find: 
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ев (2.15) 
Xo 
Similarly as in the case (2.14) for 0 «x < x, x, being defined by de MR pn 
Xo Xe 


circular motion is impossible. This fact might be connected with the repulsive forces 
for very small distances suggested by some authors (Werle 1953, Szamosi 1953, Marx 
and Szamosi 1955) attempting to explain the „Jastrow well“ (Jastrow 1951) with the 
help of „К. r.“ for the scalar field. 
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(е) The rigorous “k. г.“ effects in the scalar field theory. The results of the pre- 
ceding case (d), may be easily shown to be unvalid in our two body problem. Namely 
жое 2 
the neglect of the expressions Vi — В? in the equation of field together with x a 
с t 


is not justified. The apparently “consequent neglect of all the terms ~ as in this 
с 


equation misleads іп two different respects: (i) the neglect of the finite velocity of 
propagation of interaction, which for very small distances seems intuitively to be 
а good approximation and (ii) the neglect of some aspects of „К. г.“, the important 
features of which are the coefficients V1— f? in the sources of ф(х). 
2 

Now we shall show that “f. г.“ E a being neglected and үз — pi being 
retained in the sources of ф(х) quite different qualitative results from those of the (9) 
case will be obtained. Our equations of motion are: 
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Similarly as in the preceding cases we find for our circular motion? 
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From the character of these equations it is seen that for arbitrarily small «a circular 
motion is possible. à 

(f) Comparison of the results for all the cases. Fig. 1: В? = v?c? as function of x 
for all the discussed cases for the two bodies undergoing a circular motion on a circle 


of diameter г = ~, where и = 0 in the field equations. The curve 1, corresponds to 


the nonrelativistic case, 2, — the zero cómponent of a vector case with neglect of 
„f. r.“, 3, — the case of „К. r.“ with the additional neglect of relativistic terms in ae 
sources of the field, 4, — the case of the neglect of „f. r.“ only, 5, — the case of „К.т. 

and „f. r.“ being both retained (the rigorous solution). In Fig. 2 we have the correspond- 


E 
5 For the sake of quantitative discussion see (e.g., Figs. 1 and 2), we shall approximate = 


by (the „Hulthen expression“); this approximation is very good for very small zz, as it is 


1—ех 
the case in our future discussions. 
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Fig. 1. В? = v?/c? as function of x = pur for all the five cases discussed in Sections 2 and 3 for y = 0. 


Curves 1,: non-relativistic, 2}: vector type interaction, Зу: scalar, with neglect of Vis B;? in the sources 
of ф(х), 4: scalar, with full kr", 5,: scalar, with full „К.т.“ and „f.r.“. 


Xc 0.05 56 
Fig. 2. В? = v?/c? as function of x = иг for all the five cases discussed in Sections 2 and 3 for u Æ 0. 


Curves 1,: non-relativistic, 2,: vector type interaction, 35: scalar, with neglect of V1 — Bj? in the sources 
of ф(х), 4»: scalar, with full „k. r.“. 


ing functions f? = В?(х) for и-Е0 i п the equation of field except the curve for the 
rigorous solution (this last would be connected with great mathematical difficulties). 
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The differences between our curves are important only in the region of very small 
x (0 & x <0.1). For large x all the curves converge to the curve for the nonrelativistic 
approximation. All the relativistic curves except 3, and З, which come up to x, only, 
give В = I for x = 0. 

From Fig. 1 and 2 the differences between the respective magnitudes of В? for 
all the cases may be estimated and thus the validity of the respective approximations 
may be appreciated, 


The differences between the curves „n,“ and „A are important for x 2» x, 
where the curves „лу“ decrease as 1/x and „пз“ exponentially.— It is hoped that the curve 


for the rigorous solution foru=F0 would be similar to the curve 5,. In the region of small 
x the factor e" is unimportant. Only the characteristic not renormalizable “tails“ 
of selfaction might alter the qualitative result. However according to the result of one 
of us 6) the selfaction in the case of circular motion for the symmetrical Green's func- 
tion causes additional repulsive (excentric) forces smaller than the effective interaction 


forces (by the factor —2х, 0.028). 


3. The Rectilinear Motions and the Effective Forces 


The repulsive forces for small distances resulting from “k.r.“ for one particle 
in a given external scalar field were discussed by. Werle who suggested that 
this result remains also valid for two bodies interacting through a scalar meson 
field. It seems at first sight that this cannot be so, as the radii of the circular motions 
are arbitrary’. However, every point of the curve 5, corresponds to a different energy. 
From the fact that the circular motion with a given, arbitrarily small radius is possible 
for some definite energy it may not be concluded that for some other energies the 
forces between the particles may not become repulsive. The result mentioned in 
82 states only that for small distances the repulsive forces may not be understood 
as a rigid barrier between the particles, below which stationary motions are impossible. 

In the present section we shall investigate the effects of “k.r.“ in the case of 
rectilinear motions for some definite energy. This method will enable us to approach 
the problem of the effective forces between the particles in the scalar meson field 
theory. 

(a) Effective forces. Let the motion 7(t) be described by the differential equations 


for which we have the following conservation principles 


e = Е(у°,г) (3.1a). 
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6 J. Sawicki, Acta Phys. Ас. Sci. Hung. in print. 
7 We areindebted to Professor J. Werle for discussions on the subject. 
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йг - = 
where v = An |r|, e and J are constants. 


Let m in the equation of motion denote the rest mass of the particle. Naturally 
we shall call the „effective force“ the quantity 
Fy = m 5 for definite e and 17 
This quantity may be obtained in some different forms with the help of the inte- 
grals of energy and angular momentum. Let us suppose that (3.1а) may be solved, 
for vas v? = ylr,e). We shall denote f(r,s) = g(y,r). It may be easily shown that 
Ey is у 
тд l = 3 
== 3.2 
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The last expression is an Е of the Lorentz force, F,, may be written also 
in the form? 
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where 4 = Ie may be expressed by r with the help of the equation 
y = J2fr + Ar? 


We may also write: 


Fy=¥V (ae n (fr i$ye). т (3.2 c) 


The last equations are useful for the discussion of the relaiivistic effective forces 
between the particles. 
(b) One body in an external field. Let us consider the particle of mass m and 


position vector & (t)in the external meson field ф(х) for which Ф, = 0. According to Eq. 


тс? — 
(1.3) there exists the energy integral Ё = == . According to Marx and Szamosi 


(1954) the equation of motion may be transformed with the help of this integral into 
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the form: m = = VU, where U= lep _ (ev) | LS . U is the rigorous effec- 
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tive potential; the term — 
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is connected with the repulsive forces resulting 


from „k.r.“, if ф exceeds mc, If p = g(r) the last formula may be directly obtained 


А .mol- mol- 
* However, in general div ri — 15 0 and — — — I may not be represented by the curl of a vector. 
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- dr И ; 
8 Since А = r = 0 in the circular motions, the last form would be convenient for the 


discussion of the effective forces in these motions. 
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with the use of (3.2а). In this case we have the angular momentum integral 
т — 81сФ z Ех dé dé 
yl— в — p? dt 
of the energy integral. Since f= © =const. and 9 = 282 = с? = G 
c € 
we obtain the effective force of Marx and Szamosi with the help of (3.2). 
Let ф = g/r; let us consider the one-dimensional motion: $, == x(t), £ = £, = 0. 


With respect to the energy integral we have here: €= Le. 2° 


TE па me: 

Є becomes determined by the initial conditions. 
In Fig. 3 the function В = f(x) is given for different values of €. All the curves are 
focalized in the points where р? = 1, |x| = r,. On expressing Pin terms of x or xin terms of 


ET which may be rewritten as ғ/с E эзе J with the help 
t 


D and on integrating one obtains the time dependance t = | dt = E T However, 
we observe that these equations contain some contradictions. Namely for the initial 
conditions characterized by Є > 1, By < 0, хо > г, (the „direct fall“) our point reaches 
the point В = —1, x = г, at t = со (the respective arcs in the diagram are indicated 
with arrows). In this singular point, mathematically, the history of the motion is 
finished. Physically it means that the point reaches В = —1 for x — r, and further 
nothing may be told of its motion similarly as in the classical case of the fall of the 
Moon on the Earth. If € < 0 (e.g. By < 0, 0 < x, < r,), the body would be initially 
attracted, then repulsed, and its motion would find its end in the singular point x = г, 
В = 1 after some finite time. On admitting the theory connected with these motions, 
we would rather obtain physically singular motions in the region of the singular points 
induced by „К.г.“ and not the repulsive forces. The last difficulty may by avoided 


> D D 
by replacing in the variational principle the terms — m;c f ds; + glc f ф (©) ds; Ьу! 


— f [m;c — "i ф(ё)| ds; We obtain then the equation of motion for one body in a 


given external field 
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ds 
For ф = 8. . in the case of one-dimensional motion the energy integral results in 
"neat E : ў $ а . 1 el b Dez 
the form: € —————— . The curves В = f(x) are identical with the respective 


уту. 
curves from Fig. 3. Only the arcs beneath т„ do not correspond to negative energies: 
every arc of some positive energy is nern by the respective arc of the same 
` energy beneath r, (see Fig. 4). We may interpret the motion starting with the initial 


О 


"1 
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. 10 We are indebted for this remark to Professor Werle. 
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conditions indicated with a dot in Fig. 4 as an attraction up to x = v,, repulsion beneath 
r, and then pushing the particle up to infinity. 

The present change of the theory with the introduction of the absolute values 
into the variational principle seams, ны. to be of a rather doubtful character. 
It induces the appearance of „signum“ functions depending on ф in the sources of 


K СС А 
-1 1|8 
Fig. 3 Fig. 4 Fig. 5 


Fig. 3. x = иг as function of В = vje for the „direct fall“ of one body in an external field (u = 0 in 9(7)). 
Fig. 4. x = ur as function of В = v/c for the „direct fall“ of one body in the case of the „theory with the 
absolute values“. 
Fig. 5. x = иг as function of В = v/c for the „direct fall“ of two bodies (м = 0 in g(r)). 


the field and so gives in a certain sense a non-linear character to the field equations. 
Apparently the problem of sign is not important, since the „signum“ functions appear 
only as coefficients of the д functions and ф diverges on the world lines. However, 


me gÈ) 
SS e(£) 


The same expressions in the inertial terms of the equations of motion are renormalized 


it seems that in the expression #7 m; might not be treated as definite constants. 


1 
with the introduction of the divergent term of g(5) into m,. The consequent renormali- 
zation procedure in the equations of motion and the field equations yields complicated 
equations different in form from the conventional ones. 
(c) The direct „fall“ in the case of two bodies with the retaining of „Е.т.“ in 


the scalar field theory. From $ 2 we know, that the expression = p? must be 


consequently retained in the sources. We shall discuss the one-dimesional motion des- 
22511629 IE RUE 


1 
cribed by eq.(2.16) for u = 0: £, = x(t) = — £,, £j = $, = Ё, = &, = 0. For x (t) we get 
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d 1 — 1 2?\ dx 9 g? 
= 1 GES ыш. мш, م‎ 
dt V1 — В? (m ү В С d dt a дх r (3-3) 
where 
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== 2x SS =з 
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This motion in the approximation of »k./r.“ describes the direct „fall“ of two 
bodies of equal masses for u = 011, 
From Eq. (3.3) we deduce the energy integral 


1 1 А 
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€ 
where £ — Am . In Fig. 5 В = A(x) is shown for different г. With the 


help of simple integrations the time dependence may be effectively calculated. 

From Fig. 5 it is seen that if the initial conditions are such that f, < 0 and the 
point хо, Во lies above the critical curve x = r, (1 — f?) (being not the integral curve) 
zo being finite, x reaches 0 at В = —1. Apparently no repulsive forces appear at small 
distances. If By < 0 and the point x, f, lies below the critical curve x = г, (1 — 82)’, 
the motion will be initially slowed down, and afterwards it will go towards larger 
x's up to the critical point lying in the intersection of the integral curve with the critical 
curve x = г, (1 — f). In this point (which is reached after a finite time) the motion 
is mathematically finished. Physically this is strange, because our two bodies have 
then a finite distance and finite outwardly directed velocities. It the initial conditions 
define a point on the left branch of the critical curve, an ambiguous solution results. 
If this point lies on the right branch, nothing may be said about the motion. 

It may be easily understood what is the reason for the appearance of the critical 
points at x =r, (1 — 8%). K = x — г, (1 — fi?) ^ Е 0 is the condition expressing 
the energy integral in terms of В and at the same time it is also one of the coefficients 
of the second derivative in the initial equation of the second order. 

It may be concluded that the „К. r.“ approximation is inadequate for the problem 
of the „fall“ of two bodies for small mutual distances. It is, however, characteristic 

.that the analogous difficulty for one body in an external field induced other result. 
If the initial conditions from the region dashed in the diagram inducing unreasonable 
motions are excluded from the consideration, everything is physically consistent for 
other initial conditions and then the effective forces are attractive. 


For J — O and for different £ — E the effective forces may be obtained graphically. 


E 

11 The assumption д = 0 is of no importance for the qualitative investigation of this motion. The 
same equations result in the „theory with the absolute values in the variational principle“ under the 
assumption that the singularity in the sources is compensated by the „signum“ functign in the sources. 


470 J. Plebanski and J. Sawicki 


They appear attractive up to х = 0. It seems that if repulsive forces between two 
bodies exist in the classical scalar theory, they may not be connected with „К.т.“ only. 
It is possible that after the application of „Ёг.“, i.e. after taking into account 
retardation and selfaction, the repulsion at small distances would appear analogically 
to the case of one body in an external field. The investigation of this problem with 
the use of „К.г.“ only results in physical contradictions. It may be concluded that all 
considerations with the assumption of „k.r.“ concerning the two-body problem should 
be done with extraordinary care. All the physical conclusions from ,,k.r.“ should be 
stated exactly and not only intuitively. 
We are much indebted to Professor L. Infeld for his kind interest in this 
work and Professor J. Werle for helpful discussions. 


KPATKOE СОДЕРЖАНИЕ 


Плебанский и Савицкий, ЗАМЕТКИ OB УСЛОВНОМ ВОПРОСЕ 
ДВУХ ТЕЛ В КЛАССИЧЕСКОЙ СКАЛЯРНОЙ ТЕОРИИ МЕЗОННОГО ПОЛЯ. 


В $1 поданы некоторые общие замечания, касающиеся проблемы „Cauchy” для 
системы скалярное мезонное поле + частицы. В $ 2 подано решение вопроса двух 
тел в случае кругового движения для j, —0. Решение это было сравнено с решениями 

‚полученными для других приближений. В $ 3 поданы некоторые замечания 
о прямолинейном движении, касающиеся вопроса сил отталкивающих, высту- 
пающих на небольших расстояниях между нуклеонами в теории скалярного поля. 
Доказано, что при пренебрежении затухания из теории вытекают некоторые 
физические противоречия. 
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THELOCALIZED MOLECULAR ORBITALS AND THE STRUCTURE 
OF SOME XH, MOLECULES 


By WEODZIMIERZ Коғоѕ 
Institute of Physics, Polish Academy of Sciences, Warsaw 
(Received June 10, 1955) 


The method of equivalent orbitals in the two-centres approximation for bond orbitals 
is-used to calculate the variation with bond angle of the total energy of the NH, and PH, 
molecules. A modification of the orthogonality conditions for the orbitals enables to calcu- 
late the bond angles of the equilibrium configuration in a qualitative agreement with 
experiment. It is shown that the theory in this approximation over-emphasizes the hybri- 
dization of the lone pair orbital as compared to that of the bond orbitals, and over-empha- 
sizes the part played by the lone pair electrons in determining the molecular structure. 


I 


The problem of bond angles has been treated recently by various authors. Pople 
(1950) using the method of equivalent orbitals developed by Lennard-Jones and co- 
-workers (Lennard-Jones 1949a, b; Hall and Lennard-Jones 1950; Lennard-Jones and 
Pople 1950), and approximating the equivalent orbitals by the orbitals formed from 
the atomic orbitals on two centres, has calculated the equilibrium configuration of 
the H,O molecule. Mellish and Linnett (1954) have shown that introducing the electro- 
negativity concept to the considerations a qualitative-agreement with experiment can 
be obtained for bond angles of the XH, and XY,, molecules, where Y is a halogen atom. 
The authors agree in the conclusion that the lone pairs play a significant part in de- 
termining the bond angles and the molecular structure, and that it is not justified 
to regard the bonds as formed from the pure p-orbitals of the central atom. Conse- 
quently Pople (1950) and Mellish and Linnett (1954) regard the valencies in the mo- 
lecules with lone pairs electrons as being approximately tetrahedral, the effect of lone 
pairs being such that the bond angles are reduced below the tetrahedral value. 

However Artmann (1954a) using the semi-empirical valence-bond method, and 
assuming the bond eigenfunctions to be formed from the pure p-orbitals of the central 
atom was able to calculate the bond angles in XH, molecules in satisfactory agreement 
with experiment. Taking into account the effect of hybridization of atomic orbitals 


(Artmann 1950 b) he obtained a farther improvement of his results. 
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In the present paper we have calculated the variation with bond angle of the total 
energy of the NH, and PH, molecules using the method of localized molecular orbitals 
in the two-centres approximation for the bond orbitals used by Pople (1950) and 
Duncan & Pople (1953) to determine the equilibrium configuration of the water mole- 
cule and the hybridization in some simple molecules. 

It is shown that by introducing à modification of the orthogonality conditions 
for the localized molecular orbitals the theory can give in agreement with experiment 
the bond angle in the PH, molecule smaller than that in the NH, molecule although 
the 2s—2p level separation of the nitrogen atom is greater than that of the 3s—3p 
levels of the phosphorus atom. However, it has been found to be impossible to obtain 
in this approximation in the pure theoretical way the bond angles in quantitative agr- 
eement with experiment. А. more detailed analysis of the results shows that the theory 
over-emphasizes the hybridization of the lone pair orbital as compared to that of the bond 
orbitals and consequently over- emphasizes the part played by the lone pair electrons 
in determining the molecular structure. It seems therefore in connection with the 
results of Artmann (1954 a, b) that although the hybridization is an important factor 
in determining the molecular structure, in the simple approximation it is more justified 
to regard the bond angles in XH, molecules as deviating rather from 90° than from 
the tetrahedral angle. 


II 
We assume the molecular bonding orbitals to be of the form 


X (b) = АФ (X36) +u P (Hi; 15), (1) 
and the lone pair orbital 
x0) = Ф(Х; 1), (2) 


where Ø( X; т) (m = b,,l) are atomic orbitals of atom X (X = N,P) formed by hybridi- 
zation of ns and np functions (n = 2,3 for nitrogen and phosphorus atom respecti- 
vely). 

Explicit forms of the hybridized atomic orbitals are most simply referred to axes 
of quantization taken along their direction of maximum electron density with the 
atomic nucleus as origin. So we shall have 


D(X; b;) = cos в, D(X; ns) + sin e, D(X; np(b;)), 

D(X; 1) = cos e, D(X; ns) + sin ej D(X; пр(1)), (3) 
where Ф( Х,пз) and (X; np(b;)) or D(X; np(l)) are the atomic ns orbital and the np 
orbital pointing in the direction b; or J respectively; e, and ej are the hybridization pa- 
rameters for the bonding and lone pair orbitals respectively. Taking cos e, = 0 and 


cos & = 1 we have the atomic s?p? configuration, while cos e, = 1/۷3 and cos Vus = 0 
give rise to an atomic sp* configuration. 


Tee ee 
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Let х denote the bond angle and В the angle between the lone pair and a bond. 


They are related geometrically by 
2 
sin В = —— sin Ge 


3 2 
Following Duncan and Pople (1953) we shall assume the the orthogonality conditions 
for the molecular orbitals (1) and (2) to hold. This gives the following relations between 
the angles х and В and the hybridization parameters e, and €, 
A? (cos? e, + sin? e, cos х) + 24и (S, cos e, + Sp sin e, cos a) +u?S5=0, (4) 
A (cos є cos e, + sin єр sin €, cos В) + u (S, cos ej + Sp sin €, cos В) = 0, (5) 
where 
S — 1 Ф(Х; ns) Ф(Н; 1s) dr, 
$, = | Ф(Х; np(1)) Ф (Hay; 15) dr, (6) 
Sy = | (Hay; 15) (Ha); 15) dr, 


are the overlap integrals. The directions of the p-functions of the central atom are 
assumed to coincide with the directions of the bonds. The normalization condition 
for the bond orbitals gives the relation between A and u 

A? + 2Au(S, cos e, + S, sin &j) + u? = 1. (7) 
In the calculations of Pople (1950) for the water molecule the bond polarity para- 
meter A/u was estimated from the dipole moment of this molecule. We have also 
found it impossible to determine the ratio A/u by minimizing the total energy; how- 
ever for the PH, molecule it was also not possible to determine A/u from the dipole 
moment. In this approximation it is impossible to find such a value of A/u which would 
give the correct dipole moment of the PH, molecule. This is why in the calculations 
of the energy we have made an additional assumption that the square of the ratio 
Alu is equal to the ratio of the electronegativities of the corresponding atoms, i.e. we 
have taken A/u = 1.2 for the N-H bond orbital, and A/u = 1.0 for the Р-Н bond orbital. 
It has been proved that our general conclusions are not sensitive to the value of A/u 
assumed, within reasonable limits, in the calculations. 

In order to avoid the three centre integrals in the energy calculations an additio- 
nal approximation used also by Pople (1950) for the charge distribution in the bond 
orbital has been introduced, i.e. instead of 

x (b) x (b) = 1? Ф (Xib)  (X;b) + 24u O (Xib) Ф (Hals) 
+ и? O (Hqls) Ф (Hiis), | (8) 


we shall assume 


х (6) x (b) = (A? + AS) D (X;b) D(X;b) + (u? + Au) (Н1)Ф (Нуз), (9) 
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where 
S == 5, cos & + Sp sin e. 


The right-hand side of (9) répresents the same total charge as (8) distributed in a similar 
manner. The justification of this procedure can be found in Ruedenbergs paper 
(1951). 

The total energy of the molecule was calculated from the well known formula 
(Lennard-Jones and Pople 1950) 


Lor» Hun У»! (2(nn|G|mn) — (mn|G\nm)} + Ў; ZaZglags (10) 
n m n а+в 


where the Latin letters refer to the occupied space orbitals, and Greek letters are used 
for the nuclei. For abbreviation we have written 


(Im|Glnp) = | f ж) Xn Dn © (дат; dr; (11) 
Н = fn.) H; tn (i) а (12) 


and Н; is the operator representing the kinetic energy and the energy of interaction 
with the bare nuclei 


1 
EE EE аз) 


The energy єап be divided into the following parts 
| E = Bo + By + Ey + Е + Ep (14) 


where Еъ,ь, represents the energy of interaction between the electrons and the protons 
of different bonds, Е», is the energy of interaction between electrons and nuclei the 
same bond, E, is the energy of interaction between the electrons in the lone pair 
and the electrons and protons in the bonds, E; represents the repulsion between the 
electrons in the lone pair, and E, is the kinetic energy of the electrons and their poten- 
tial energy due to the nuclei (and inner shell electrons) near which the electrons are 
concentrated. | 
Explicit forms of the various parts of the interaction energy are 


3 
Esp; = SE 12(2# + AuS) Win (n.&,*) — 12(u? + AuS)W iw (онн) 


+ 19005 + диз). (se eps) + 24(42 + ди) (и? + Aus) us. 
X Wig (п,єь,ә) + 12(u? + AuS)? ЁЎин(онн), 
15 | 
Eu = TA 6(4? + A uS) Wan (п,вь.0) — 30(u? + AuS)W un (онх) 
=: + 3(42 + Аб) „(пьвъьвь,0) + 6(22 + AuS) (и? + AuS) (16) i 
X И’ н(пьгь,0) + 3(u? + AnS)? нн (О), | 
Вы = —6Wiy(n;ey В) + 12(42 + AUS) № neen f) + 12(и? + А5), inn, (17) | 
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Ey = W,2(n,€;,€),0), (18) 
E, = 2cos*e,W,, + 2 sine W np + 6(A? + AS) cos?e,W,, 
+ 6(4* + Au S) sin?e,W,, + 6(u? + AUS) Wap (19) 
where 
2 
Wy 9(1,8},82% з) = De Ф1(1) Ф; (2)(1/n ) тат» (20) 


represents the interaction between two electrons in the atomic orbitals Ф, and Ф, 
formed Ьу hybridization of ns and np atomic functions with hybridization parameters 
е; and =», on the same centre and inclined at an angle a, ;; 


Wr) = || 910) 0232) Und)dudr, (21) 


gives the interaction between electrons in a hydrogen 15 orbital and a hybridized 
atomic ns—np orbital on another centre, whose direction is inclined at an angle ox 
to the line joining the centres; in 


lp) = f 410) (Шт) (22) 
the electron in the hydrogen 15 orbital from (21) is replaced by a hydrogen nucleus; 
Wann) = | f $1.) 2,02) (1m ddr, (23) 


represents the interaction between two electrons in the Le hydrogen orbitals separated 
by the distance онн; and й 


W nun) = T pall) (ln ld GA 


gives the interaction between an electron in the 15 hydrogen orbital with a hydrogen 
nucleus at distance Oyy- 

° The terms W,, and W,, represent the kinetic energy of ап eletron in the atomic 
ns or np orbital respectively and its interaction with the nucleus of the central atom 
and the inner shell electrons; №, represents the kinetic energy and the energy of 
interaction with the proton of an electron in the atomic Is orbital of the hydrogen 


atom. 


In the terms representing the interaction between the IST or the electrons 
in the 1s hydrogen orbitals with the inner shell electrons the latter have been considered 
as point charges situated in the nucleus. 

Following Pople (1950) we have neglected the exchange energy terms. aE 
as may be seen from (13), we have taken into account the „local ery given by 
7 Е, + Е, + Ey which was neglected by Pople. 

` The integrals (20) to (22) are evaluated in the Appendix using Бре (1930) 
approximation for the atomic wave functions of the central atom. Only for Wa, and 
Wop in the case of NH, molecule we employed more exact values obtained by Moffitt 
(1950) using the functions of Duncanson and Coulson (1941). The integrals (23), 
(24) and (20), (21), (22) for n = 2 were evaluated by Pople (1950). The distance 
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Х-Н was assumed to have a constant value 1.01 A and 1.42 A for the NH, and PH, 
molecule respectively, which are the experimental values for the equilibrium configur- 


ations of the molecules. 


III 


The results of the calculations may be summarized as follows: 

(1) Assuming both orthogonality conditions (4) and (5) to hold, the energy has 
been found to have no minimum. 

(2) Assuming instead of (4) and (5) the orthogonality conditions for atomic orbi- 
tals of the central atom: 


cos?e, + sine, cos х = 0, (25) 
cos £j COS €, + sin е, sin €, cos В = 0, (26) 


which are obtained by putting u = 0 in (4) and (5), a deep minimum of the energy 
has been found near a = 102° and a = 98° for the NH, and PH, molecule respectively. 

(3) The orthogonality conditions (25) and (5) were used by Pople (1950) in his 
calculations of the equilibrium configuration of the H,O molecule. In this approxima- 
tion we obtained a minimum of the energy near « = 110° for the NH, molecule and 
х = 106° for the PH, molecule. The experimental values of the bond angles are 
a = 106?47 (Herzberg 1945, p. 439), and х = 93°50’ (McConaghie and Nielsen 1953) 
for the NH, and PH, molecules respectively. The variation with bond angle of the 
various parts of the energy for the NH, molecule is shown in Table I. The hybridiza- 
tion parameters of the bond and lone pair orbitals are also given in Table I. 

(4) The failure of the calculations reported in (1) is due partly to an over-emphasi- 
zing in the orthogonality conditions of the overlap between the hydrogen 1 5 functions. 
This may be seen from the fact that the minimum value of the bond angle corresponding 
to cos £j = 0 for the PH, molecule, obtained from (4), is « = 95°20, i.e. greater than 
the experimental value. However Coulson and Duncanson (1943) have shown that 
a larger value of the effective nuclear charge must be taken in the atomic orbitals 


TABLE I 


104° 


120° 


cose; 0.19480 0.25486 0.26383 0.33333 
сов?є 0.24634 0.14170 0.12589 0.00000 
Eu 1.9952 1.9738 1.9697 1.9271 
ER 0.9021 0.8879 0.8868 0.8803 
Ej 3.9686 4.0031 4.0127 4.1243 
Ey 0.8558 0.8250 0.8192 0.7622 
E, —17.7757 — 17.7602 —17.7585 —17.7401 
E — 10.0540 — 10.0704 — 10.0701 —10.045 2 


` : 
Әһ 
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which form an attractive molecular orbital than in the orbitals of the isolated atoms. 
The influence of the variation of the hydrogen »screening constant“ Z н on the dipole 
moment of the ammonia molecule has been also considered by Duncan and Pople 
(1953). So we have repeated the calculations varying 2н in the hydrogen 15 functions 
in the orthogonality conditions (4) and (5). 

In the case of the ammonia molecule the results are 


Zw 1.6 
TER 


1.95 2.2 3.3 


1102 108.5281057 


For the phosphine molecule we obtained 
Zr AP 3.5 5.0 2276-5 


eo. МБ МОУ WS 


The screening constant giving for the NH, molecule the bond angle in agreement 
with experiment has the value Zu = 2.5 and for the PH, molecule Zu = 5.3. 

The variation of various parts of the energy for the ammonia and phosphine 
molecules calculated assuming Zy = 2.2 is shown in Tables II and III respectively. 
The significant part played by the orthogonality conditions for the orbitals in determin- 
ing the minimum can be seen from the results. 


TABLE II 
0.008447 0.14214 0.21462 0.26841 0.31554 
cose, 0.81320 0.40434 0.23161 0.10800 0.00000 
Ey, 2.0165 2.0175 1.9971 1.9730 1.9461 
Lë 1.1401 0.9672 0.9462 0.9404 0.9395 
Es 4.2328 4.0018 4.0064 4.0586 4.137 
Ey 0.8054 0.8787 0.8522 0.8124 0.7632 
E, —18.1131 —17.9097 — 17.8578 —17.8245 ° —17.7972 
Е — 9.9183 —10.0445 —10.0559 —10.0401 —10.0112 
TABLE Ш 
92° 96° | me 110° 120° 
cose, 0.008748 0.05534 0.16493 0.23944 0.34116 
соз?ву 0.84116 0.67855 0.40303 0.23226 0.00000 
Dë? 1.1381 1.1408 1.1310 1.1133 1.0765 
£s 0.1150 0.0587 0.0185 0.0100 0.0096 
Ej; 2.3780 2.3108 2.2510 2.2688 2.3655 
Ey 0.4655 0.5008 0.5183 0.5024 0.4479 
E, —12.4297 ; —12.3667 — 12.2907 —12.2503 — 12.1960 
E — 8.3331 — 8.3556 — 8.3719 — 8.3558 — 8.2965 
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A detailed analysis of the dependence of the various parts of the energy and 
bond angles on the hybridization parameters shows that the theory over-emphasizes 
the hybridization of the lone pair orbital as compared to that of the bond orbitals and 
consequently the theory of equivalent orbitals in this approximation over-emphasizes 
the part played by the lone pair in determining the molecular structure. 


Appendix 1 
The integrals (20), (21) and (22) for n — 2 were evaluated by Pople (1950). For 


an arbitrary п using Slater's (1930) approximation for atomic orbitals 


Zxr 
D(X; ns) = Ме 


e Zur 
D(X; пр) = №, V3co$ де, (А,1) 
where 
2 Д\?п+1 1 
Nes E | Art 2n!’ (4:3) 
we obtain 


Win (пам) = N? EN + 23 sin Е cos € COS 01у A, (Z,0,n,1) 


(A, 3) 
+ 3 sin? e P,(cos ay) E (Z,0,n,1) — > A, (2.0.5 alt, 
where 
2(n—1) 2Zr, 
q n—1) e-—— — mr, 
A,(Zri,n,p) = [| сова er” de: (А,4) 
ry? 
N e | 
7 н(пь=, ан) = m {Ton(Z) + 2V3 cos e sin e cos “yg In (Z) (А,5) 
i+ 3 sin? e P, (cos ан) Го» (Z)}, 
— = IZ) = ER (72,0, nl) — A,(Z,2,n,1) — A,(Z,2,n,0), 
Е | 
7T fes (Z) E A, (Z,0,n,1) Te A,(Z,2,n,1) ZS ER (Z,2,n,0), (A,6) 
1 I 3 р 1 
g 2" (2) = CH [A ,(Z,0,n,1) — As(2,2,n,1) — A,(Z,2,n,0)] — эл Тол (2); 
Маа)? - | 
-Wia(n, 8,835,045) = №(4л)? 1B(1,n.Z) + — cos & cos ез sin €, (A,7) 


X sin =, cos & 5 B(2,n,Z) + 25 ain & sin? £5 Р, (cos х) B(3,n,Z), 
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where 


m—1 
B(m,n, Vr) == [fr г" е = ( li FE To ) E dt, dt, (A,8) 


T> 


The explicit expressions for the integrals A „(Zm;n,p) for п = 2 can be found 
in Coulson’s (1942) paper. The integrals for n = 3 are given in Appendix 2. 
For n = 2 we shall have from the formulas given above! 


’м(2, вм) = —— 2 МЕ, 0,0) Е Ј (2,0,0) cos ғ sin € cos ам (A,3) 


Aas 


25 D 
-Е 32x E (Z,0,0) EX E i sin? e P, (cos ci 


Lr А 
W н(2,8,01н) = 067i Го,» (Z) + —— Aa (2) cos e sin € cos ayy (А,5') 


= 2 


І, (Z) sin? e P, (cos xyz), 


75 
+ 48л? 


where 


=: h,a (2) = 5(2,0,0) — (2,2,о) — Ль(2,2,0), 


1,2 (2) Fr. Ja(Z;0;0) x Ja (2,2,0) ‘>, Jao(Z.2.0). | (А,6’) 


3 1 
L 2 2( =È Un(Z09) — Jay(Z2,0) — Ia(2249] — 5- 1400), 


93 185 
EOE peto з) = 2515 du 1159 sin & sin Ca COS € COS £9 COS 05 (АТ) 
E Ar sin? e, sin? e Р, (cos gl 
and for n — 3 | 
W,nß,&&n) = gi Ja 0,0) N es 51 Ј,„(0,0,0) cos e sin e (A,3”) 


1 
X cos ум + Lus sin? e P, (cos an) KG 0,0) — E: Legal 


where 


1 We give here also the formulas for п = 2 because there are some misprints in the formulas 
У 
published by Pople. 
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ЛО ا‎ == == 


w’ 2 
Wig (3,8%) = 173-6 Го, (6) + E 1,3 cos € sin € COS ur (A,5”) 
©? P 
+ 333-61 I, del sin? e P, (cos ин) 
where 

Ja. ol = Лаа(0,0,0) — J44 (0.2.0) — Jas (0,2,0), 
d ; 
pu T, del = Ja, (0.0.0) ж J47(@,2,0) be Jas(%,2,0), (A,8’’) 


| 3 1 
л Г,.з(о) = 9 [J49(0.0.0) SS Jag(0.2,0) = Jalo , 2.0)] "E Эл То, sto, 


793 2979077 , 
Я. ==> = ZI —— 27 sin & sin €, COS Eu COS £9 COS @ 
12(3,61 s€ 25% 2) toate 34560 1 2 1 2 12 (A.7”) 
2 
SE a sin? e, sin? ee P, (cos оз) - 


Appendix 2 


The integrals given below, evaluated by standard methods, may be regarded as 
a continuation of those published by Coulson (1942). The notation is the same as 


in Coulson's paper. The convenient abbreviation 


X= «? — ур 


is also used. 
a 
Je, Bäi = fe e e Pre dt 
D 
= E (2402 (32^ + 10426? + 385)e-/e — [05X5 + 10* Eis 
+ 2092 X*(3a? + В?) + 2409? X?a(a? + В?) + 12005 
+ 100262 + 8%) + 240a(3a* + 10028? + 38*)]e-*e), 
AE " rers P dr 
8 | 
= 25, (о Х(8о4 + 1028? + 38%) — 88(To* + 14:98 + 389] 


x 120ae-fe + [05X5 + 2004 Хо + 4003ХЗ(5а? + В?) 
+ 24002 Ха (ба? + 382) + 1200Х(35а4 + 4202 6? + 308%) 
+ 960a (Tat + 140282 + 38*)] Bee}, 
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DECHE f: Ta COS 0, е-“та—Вть {т 


= 2 {24[0° Х(35а* + 420282 + 384) — 24(о8 + 1) (2104 


+ 18a*8? + В*)]е-ве + [96 X* + 2095 Xo, + 120*X3(17&2 + 3?) 
+ 480% X*a(27a? + 138?) + 24 Xo*(217a4 + 174028? + 954) 
+ 576(0 + 1) (2104 + 188282 + 85)]e-«ey, 


4 
S 
(ро) = {pe cos 0, e “ars dr 


Cle —ÀMÀ —— 1. 


= e (8ай + 42022 + 389 (До + 1) 48е-#е — am 
0 


+ 10g* Хо + 201 X*(31a? + 982) + 169? X?« (17a? + 138) 
+ 24 (35a* + 42028? + 384) (9? X + 20% + 2)]e e}, 


Jas(o, B,0) = (us cos? 0, e a Вт dt 


= SC (24« [08 X(5a + 540282 + 2184) — 8928(7a4 + 860282 
4 سے‎ 2 


+ 276“) — 480ß (Be + 1) Gei + 8?)] e fe + (97 X5 + 2008 Хо 

+ 169? X*(13a? + 26%) + 489*X?« (29а? + 118) + 2492 X 

x (261a4 + 134028? + 584) + 192о®« (07a* + 26220 — 381) 

+ 11520% (хо + 1) (Ba? + 8%] Be, 
p ве = Ji "а cos? 0, еа В dr 


x дуз (leet Stage 5 219 w 2 = T dei Lä al 


x mo cum a +5 34) Бузи | » Ji 
OAM istos, 2: ei 
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Видоизменение условий ортогональности для орбит дает возможность вычи- 
слить валентные углы конфигурации равновесия. Получается качественная согла- 
сованность с опытами. Доказано, что теория в этом приближении дает чересчур 
большие величины на гибридизацию орбит свободной пары электронов в срав- 
нении C гибридизацией орбит связывающих и приписывает слишком большую = 
роль свободной паре электронов при определении структуры молекул. 
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NEW LOW DEAD CAPACITY CONDENSERS FOR MEASURING 
DIELECTRIC CONSTANTS OF LIQUIDS 


By M. Jezewsx1 
Institute of Technical Physics, Mining Academy; Krakéw 
(Received July 3, 1955) 


A new construction of cylindrical condensers for measuring dielectric constants 
of liquids has been described. Several condensers of different capacities up to 140 pF 
have been constructed. An absolute measurement of the dielectric constant of commer- 
cial, chemically pure benzene has been made, and then by means of this benzene the 
dead capacities of the above mentioned condensers have been determined. It has been 
proved that they had a dead capacity of the order of 0,1 pF and extremely small losses. 
The volume of the liquid necessary for the filling of the condenser amounted їо 10—15 cm2. 


Although the methods of capacity measurement have recently undergone consi- 
derable improvements yet the modes of construction of the test condensers for liquids 
as well as for solids are still rather underdeveloped. The chief difficulty consists in 
the fact that for making an absolute measurement of the dielectric constant a condenser 
is necessary in which the outer coating surrounds entirely the inner one, and during 
each measurement the whole space where the electric field is being established must 
be completely filled with the examined substance. This is prevented by the necessity 
of connecting the condenser coatings to the electrodes of the test circuit. Therefore 
a certain part of the condenser is never filled with the dielectric and constitutes the 
so called dead capacity. In older type condensers the value this capacity amounted to 
several percentages, sometimes more than ten, of the whole capacity. In the Nernst type 
condensers this difficulty was overcome by a double measurement at two different 
capacities of the condenser. For different reasons this type of condenser is out of 
use; recently cylindrical condensers are most frequently used. From time to time there 
appear in the literature descriptions of new constructions of condensers for examina- 
tion of the dielectric constant of liquids (Hartschorn and Rushton 1939, Wood 1952, 
Baxter and Vodden 1952). The authors try usually to reduce the dead capacity as 
much as possible. 

The dead capacity includes first of all the capacity of the wires connecting the 
condenser with the test circuit. This capacity was being tried to be eliminated or 
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determined by different means (Hartschorn and Oliver 1929); it is of the order of 
at least several pF. As the measuring condensers have usually a capacity of the 
order of 100 pF, when filled with air, the capacity of the connecting wires cannot be 
neglected. 

The present paper describes an attempt to construct a new type of condensers 
for measuring dielectric constants of liquids, their dead capacity being as low as 
possible. Measurements which permitted to determine the magnitude of this capacity 
are given below. 

The construction of the condensers is shown in Fig. 1. It is practically a cylindrical 
condenser. Its upper part g is made of Jena glass. An empty metal cylinder d fixed 
on a glass tube f of 10 mm diameter by means of araldit! constitutes the inner coating. 
This cylinder consists of three parts put together and soldered. A metal cylinder c 
cemented to the outer glass tube constitutes the outer coating. Unfortunately it was 
not possible to fix it also with araldit, as larger glass tubes, connected with metal ones 
by means of araldit, burst. Therefore a cement of lead dioxide with glicerine has been 
used. In order to make sure whether this cement would not dissolve in the examined 
liquid, its solubility in various solvents, namely: benzene, toluene, ethyl alcohol, 
methyl alcohol, acetone, ethyl ether, carbon tetrachloride, chloroform, machine 
oil and dioxan has been investigated. It was found that the above mentioned cement 
does not dissolve in any of these solvents in perceptible amounts. 

Tube e serves to introduce the examined liquid into the condenser. 

During the measurements the outer coating of the condenser is earthed. Pin bb 
serves to connect the inner coating to the test circuit. In the upper part of the metal 
socket k there is a stopper m. By raising the pin bb and rotating it by 90° the pin can 
be disconnected from the inner coating. By an adequate rotation the stopper m is 
pushed down the support p, the pin is moved down in consequence of the pressure 
of the spring r, and its pointed end comes in contact with the inner coating of the 
condenser. The glass tube a guides the pin, preventing any side movements which 
could change its capacity with respect to the surrounding. The whole upper glass part 
of the condenser is covered with a metal net, due to which the pin capacity with re- 
spect to the earth is fixed and remains invariable. The upper part of the pin b moves 
inside the metal tube s, owing to which the movements of the pin do not influence 
its capacity with respect to the earth. For measuring purposes the condenser is 
filled so that the whole inner coating d is immersed in the examined liquid. The 
liquid does not cover only the immersed part # of the inner coating. It was proved 
by measurements that losses of the described condensers are imperceptible. 

It could be expected a priori that in a condenser of the above type the dead capacity 
would be very small, but this had to be proved. There are two possibilities of doing 
it: first, an indirect one, by measuring the dielectric constant of thesame liquid by means 


of condensers of different capacities; in the case of the dead capacity being negligible, 
WT ME EA Lr ne E Е c 


1 А cement for binding metals, product of CIBA Aktiengesellschaft, Basel. 
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the results should differ very slightly. The second possibility consists in the measure- 
ment of the dielectric constant of a liquid of known specific inductive capacity. From 
the difference between the actual dielectric constant and that determined by measuring 


Fig. 1 Fig. 2 


Fig. 1. c — outer coating of the condenser: d — inner coating; g — and f — glass part of the condenser: 
b — switching pin; a — guide of the pin; m — stopper on the pin; р — support; г — spring pressing 
the pin; e — tube for introducing the liquid. 

Fig. 2. a — outer coating of the condenser, a double cylinder; b — inner cylindrical coating; c — outer 
cylinder of the middle coating; d — tube for pouring the liquid. 


it with the help of the new condensers it is obviously possible to calculate the dead 
capacity. 

Both methods have been used. Six condensers have been built. The smallest one 
had a capacity of 10 pF, the next ones had capacities of 14 pF, 42 pF, 59 pF, 105 pF, 
140 pF. As a control liquid commercial, chemically pure benzene was chosen, which 
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shows extremely small losses. An accurate bridge method has been elaborated, by 
means of which the dielectric constant of benzene has been measured in different 
condensers. The differences proved to be very small. However, a systematical decrease 
of the obtained results for C/C, (C — capacity of condenser with benzene, C, — with 
air) with increasing condenser capacity could be seen. Hence it could be concluded that 
there exists a perceptible dead capacity. In order to determine it, it was necessary 
to measure the capacity of the condensers filled with air, and then filled with benzene 
of known dielectric constant. Let the capacity of the condenser filled with air be Co, 
the capacity of the same condenser with benzene—C, and the dead capacity — Cj. 
Then from the formula 

C— Са _ 

BEN, 
it is possible to calculate С, A difficulty lies in the fact that the dielectric 
constant of benzene, particularly of commercial benzene, is not exactly known. 
Figures given in literature for specially purified and dried benzene differ 
substantially, more than those obtained in my measurements carried out with the help 
of condensers of various capacities. In addition it is not known how far small amounts 
of impurities of commercial benzene affect the magnitude of the dielectric constant. 
There was therefore no other means than to make absolute measurements of the die- 
lectric constant of benzene. To prevent any change in its humidity during measure- 
ment, I used benzene which was standing for a long time in a large bottle, partly filled 
with air and frequently opened. For absolute measurements a cylindrical condenser 
shown in Fig. 2 was constructed. Its outer coating consisted of the double cylinder a, 
which was earthed during the measurement. The cylinder bb, situated between these 
cylinders could be pushed in and out from cylinder c and in this manner the capacity 
of the condenser could be changed. The measurement consisted in: (1) measuring the 
capacity С, of the condenser filled with air, (with cylinder bb pushed up); (2) measuring 
the capacity С, of the same condenser with benzene; (3) measuring the capacity of 
the condenser with a pushed down cylinder, the condenser being filled with benzene; 
(4) measuring the capacity C, of the same condenser after its filling with air. The 
measurements were performed in the above given order. The dielectric constant 
was calculated by the formula 


PRU: 


€ 


Аз previously mentioned measurements showed that the dead capacity was small, 
in order to determine it, it was necessary to know the dielectric constant of benzene 
with a great accuracy, at least to the third figure after the point. In order to attain 
so high an accuracy it was necessary to maintain the free surface of benzene on the 
same level during both measurements (2) and (3). This task was made difficult by the 
evaporation of benzene, as it was not possible to close the condenser hermetically. 


Вла: 
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In order to keep this level at a constant height, the following device has been used: a tube 
d, Fig. 3 was constantly supplied with benzene flowing through a capillary tube f from 
a small container, while lever g automatically maintained the level of benzene ata cons- 
tant height. On taking still some more precautions, the obtained results differed 
only in the fourth figures after the point. For the dielectric constant of chemically 
pure, commercial benzene, neither purified nor dried, the following result at the 
temperature of 18°С (with respect to air) has been obtained 


ё = 2,2866 + 0,0007. 
Then by means of the formula 


2 бб 


(where C, is the capacity of the condenser with air, C the capacity of the condenser 
filled with benzene, and C, the dead capacity) the dead capacity of condensers No 3, 
4, 5, 6 was calculated. (The capacities of condensers No 1 
and No 2 were too small, and consequently the accuracy of 
the measurement too low). Results for various condensers 
differed somewhat, but they all varied within the limits 0—0.2 
pF. For condenser No 5, the dead capacity was practically 
equal to zero, as the dielectric constant calculated with the 
help of the simple formula e = C/C, was quite compatible 
with the result obtained with the use of the above described 
cylindrical condenser. For other condensers the dead cap- 
acity proved to be negative. This can be explained by the 
fact that between the end of the switching pin b and the 
inner coating there exists a certain insignificant capacity, Fig. 3. e o СҮ 
which disappears at switching in the inner coating. up d e 
As it is not difficult to build a condenser of a capacity hich benzene flows ut 
of 100 pF and more, it is possible not to take into account 5. Jever through which 
the dead capacity in condensers of the described type, so ^ excess of liquid pours out. 
long as measurements with an accuracy not exceeding 0.2% . 
are concerned. If a higher accuracy is required, however, the dead capacity must be 
measured. А complete elimination of the dead capacity is impossible and will 


probably remain so for ever. 


KPATKOE СОДЕРЖАНИЕ 
Ежевский, НОВАЯ КОНСТРУКЦИЯ КОНДЕНСАТОРОВ ДЛЯ ИЗМЕРЕНИЯ 
ДИЭЛЕКТРИЧЕСКОЙ ПРОНИЦАЕМОСТИ ЖИДКОСТЕЙ. 
Описана новая конструкция цилиндрических конденсаторов для измерения 


диэлектрической проницаемости жидкостей. Изготовлено несколько конденса- 
торов различной емкости до 140 пикофарад. Измерена диэлектрическая прони- 
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цаемость технического, химически чистого бензола и при помощи этого бензола 
измерена мертвая емкость. Оказалась она незначительной, 
рады, диэлектрические потеры ничтожно малы. Объем ж 
конденсатора составляет or 10 до 15 смз. ў 


порядка 0,1 пикофа- Е 
идкости для наполнения 
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LETTERS TO THE EDITOR 


THE RETARDING POTENTIAL AS A STABILIZING FACTOR FOR THE 
WORKING CONDITIONS OF ALUMINIUM G—M COUNTERS 


В. SUJAK 
Institute of Experimental Physics, Bolestaw Bierut University, Wroclaw. 
(Received June 14, 1955) 


It is well known that aluminium walled G—M counters show a marked instability 
in working plateau and background, during fast rate counting ушт. 1949, Loeb 
1950, p. 482). 

In these works on the excited emissivity problem, it is proved that the instability 
phenomenon, orten called “counter self excitement“ is due to the excitement of the 
cathode (Kramer 1951, 1952; Sujak 1953, 1955; Bohun 1954, Lepper 1955). Counter 
cathode exciting seems to be looked on as a production process of some donor centres 
with a rather low effective work-function in the oxidized metal layer or in the layer of 
the decomposed quenching component of the counter filling (Sherb 1949, Nassenstein 
1954). The nature of the centres seems to be like the F- or F’- centres in alkali halides 
(Sujak 1953, 1955; Bohun 1954). 

The excited cathode of the counter can be a source of slow exo-, photo-, or thermo- 
-electrons during and after fast rate counting or continuous discharge (Kramer 1951, 
1952; Grotowski, Hrynkiewicz, Niewodniczaüski 1953; Sujak 1953, 1955; Lepper 
1955). The author represents the point of view that this is the main reason for the 
large number of spurious counts, and strong influencing of the working conditions of 
the aluminium counters. 

Regardless of the marked self excitability mentioned, the aluminium walled 
counters are very often used as thin walled ß-detectors. It became evident to the author 
after an analysis of the counter-self-excitement problem that when the counter cylinders 
are made of easily oxidized metals, grids with a negative potential against the cylinders 
should stop the undesired influence of excited cathodes. Some experiments were 
carried out with aluminium walled counter including a small mashed brass grid, 
such as used by Dyachenko (1940), to test this point of view. In these experiments the 
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grid counter was operated in the same way as that used in radiocarbon technique 
(Libby and Lee 1934, 1939). The potential of the earth was put on the grid whereas 
the cylinder could have a negative or positive potential against it. 

The construction details of the grid counter used are evident from Fig. 1. 


brass grid Al cylinder 


Plexiglas 


DEER өө з» 


nylon sphere 
= 


Tee ove 


ооо ооо о 5 


Fig. 1. The schematic view of an experimental grid counter. 


To get rather bad working conditions for the experimental counters the aluminium 
cylinders used were especially oxidized in air for 3 hours at about 500°C. Furthermore 
the counter filling contained about 20% of air besides the mixture argon alkohol: 


7 and 1 cm Hg. respectively. 


| 
9-0۷ | ху 
grid=-105 V 


0 
1020 1050 1080 1100 1150 V 


Fig. 2. Examples of the gridded counter 
characteristics when grid on different retar- 
ding potentials. 


Sec di Operating voltage 
5205 = 1080 V 
ds 
SE 
250 Sé —ł "Es 
кере к= - 200 
150 
100 


retarding potential | drawing potential 


100 60 AN 240 60 100 Volts 


Fig. 3. Dependence of counting rate at constante у-гауз 
source and counter operating voltage, as a function 
of a grid retarding and drawing potential. 


On Fig. 2. are given some examples of working characteristics of the experimental 
aluminium counters used when operating in the new manner. A strongly stabilizing 
influence on the counter working plateau can be seen, having a negative potential of 


the grid against the cylinder. 


The experimental counters have 2 cm effective length and 1 cm diameter (grid). ` 


The internal diameter of the aluminium cylinders used is 1.4 cm The experimental 
counters showed about 12 counts per minute of background, when operating with 
a 105 V negative potential against the aluminium cylinder. 


Stabilizing Factor for Aluminium G-M Counters 49] 


This manner of counter construction and operating seems to be especially impor- 
tant for the B-particle detection technique when large intensities are to be measured 
and the highest vacuum technique in preparing counters is not used. 


The author wishes to express his thanks to Professor B. Makiej and Professor 
J. Wesotowski for fruitful discusion. 


Added in proof.: After this paper was given to press A. Bohun [Czech. 
J. Phys., 5, 429 (1955)] and Г. Grunberg and К.Н. В. Wright [Proc. Roy. 
Soc., A, 232, 403 (1955)] published works on the emission of électrons due to 


the color centres in Al,0, crystals and ов the surface of abraded aluminium. 
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ENERGY INTEGRALS FOR 3D ELECTRONS 
IN BODY-CENTERED IRON 


By MACIEJ SUFFCZYNSKI 


Institute of Theoretical Physics, University of Warsaw 
(Received June 27, 1955) 


To carry out quantitatively the programme of the tight binding approximation 
for solids one requires the knowledge of the numerical values of the energy integrals. 
Fletcher and Wohlfarth (1951) in their pioneer work have calculated these integrals 
for Ni as the two-center integrals between the atomic 3d functions. The atomic functions 
have been calculated by Hartree (1936) for Cu* in the Hartree-Fock approximation. 
The example of Ni is a particularly good one as in the two-center approxi- 
mation reasonable values of the energy integrals can be obtained in such cases when 
the atomic functions of neighbouring atoms overlap only slightly; moreover, as it is 
well known, the ratio of the interatomic distance to the radius of the 3d function 
is particularly large for Ni. 

Now it was thought worth while to calculate the same energy integrals for other 
transition metals, firstly for iron. Following Fletcher (1952) analytical expressions have 
been found for the 3d atomic wave functions of the transition elements for which the 
following self-consistent calculations have been published: 
for Cr** calculated without exchange by Mooney (1939): 


RA = e е 5.46 ВЕ y | (1) 
for Mn** calculated with exchange by Hartree (1955): 
Pils (4.64 e>" + 60 e), (2) 
for Fe calculated without exchange by Manning and Goldberg (1938): 
Pa (3086 есе 50,89 2). (3) 


The analytical expressions have been found by a procedure similar to that devised 
by Löwdin (1953). They fit the numerical values very well, the discrepancy being 


6 x 107? at most. 
The similarities of:3d wave functions of different ions of transition elements may 


seem spurious if taken uncritically. Nevertheless it is very interesting to see that the 
(495) 
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exponents in the above analytical expressions are the same as those used by Fletcher 
in fitting the Cu* figures. Only the Mn** function has one differing exponent; this 
function calculated with exchange is much more compact than the others. 

The energy integrals have been calculated with exactly the same approximations 
as used by Fletcher (1952). Every integral was calculated as the difference of an integral 
over all space and an integral over the atomic sphere of the central atom. The radius 
of this sphere was always taken as half the distance between the central atom and the 
neighbour. Thisisan approximation dictated by convenience. Thus inthe body-centered 


2 ; : 1 А 
structure of the lattice constant a this radius was т V3 a for the nearest neighbours 


1 | 
апа уа for the next nearest neighbours. 


The integrals have been calculated between functions with о, л, 6 symmetry 
since Slater and Koster (1954) have given very convenient tables for reducing the 
energy integrals of the tight binding approximation to the two-center integrals between 
functions of o, л, д symmetry with respect to the line joining the central atom with his 
neighbour. 

For the body-centered iron the lattice constant was taken to be a = 2.8606 A. 
(Bozorth 1953 p. 864). The analytical expression for the atomic potential in iron was 
taken to be 


U(r) = Z,(r)/r = —(1 + 25 є-3”)/. (4) 


It must be said that this expression fits the Z,(r) function given by Manning (1943) 
for the valency electrons of metallic iron only very poorly. 

The integrals obtained for body-centered iron using (3) and (4) are: 
for the nearest neighbours 


(ddo), = — 0.69876 eV, 
(ddr), = 0.4287 eV, 
(446), = — 0.0711 eV, 


for the next nearest neighbours 
(dde), = — 0.3222. eV, 


(Фал) = 0.1593 eV, 
(ddd), = — 0.02206 eV. 


It is believed now that the absolute values of these integrals are much too great. ` 


This is because the atomic field for Fe has been calculated without exchange and it 
is known that the exchange has very serious effect on the 3d functions in transition 
elements. The 3d function of iron, especially of metallic iron (Callaway 1955), is 


much more compact than the function calculated without exchange. The above calcu- 


lation gives the evidence that even in the integrated expressions the neglect of exchange 


is not allowed. The computation of wave functions with exchange for the transition 


metals is therefore of prime importance. 
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Det e 


I thank Dr G. C. Fletcher for his kind letters, Professor J. Callaway for sending 
hic unpublished paper and Professor P*O. Löwdin for sending his papers on atomic 
self-consistent functions. | 

The students J. Baran and J. Fiutak helped me fitting the analytical expressions 
to numerical functions. 
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NOTE ON THE THEORY OF THE POLARIZATION 
OF PHOTOLUMINESCENCE OF SOLUTIONS 


A. JABLONSKI 
Physics Department, Nicholas Copernicus University, Тогай 


(Received August 11, 1955) 
A theory of the fundamental polarization of photoluminescence of solutions was 


put forward long ago (Jablonski, 1935). Recently a new form was given to the expres- 
sions resulting from that theory (Jablonski 1950), viz. 


== =1 
0 = 1k=1 | (1) 


j=l k=1 
3 3 
33) I; G;— > 1, Cr 
Bea e ZE (2) 
бл. = 3,43 , 
MES I Cu» I 
i=1 y=lk>1 


where o, and P, denote the fundamental depolarization! and fundamental polarization? 
respectively, 71, Г, and Г, are the principal transition probabilities of the absorption 
oscillator, and G,, С, and Сз the „components“ of the transition probabilities of the 
emission oscillator along the directions of the principal axes ‘of the absorption 
oscillator (i.e. along the directions of Г, Г, and Г.). 

The fundamental depolarization (polarization) is called „principal“ depolarization 
(polarization) in cases in which the same oscillator is responsible for both absorption 
and emission, і. e. when С, = Г, 6, = I, and Сз = Гу. The principal depolarization 


15 


1 The limiting value, obtained by eliminating all depolarizing factors, of the rate of depolarization 
of Juminescence emitted by an isotropic solution at right angles to the electric vector of the plane polarized 
Primary light. l А 

2 The limiting value of the rate of polarization, obtained in the conditions mentioned in the foregoing 
footnote. 
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3 2 3 З 
(xn) - п 


j=l 


[i 


0p — 3 2 3 , (1a) 
2 
| > n) +2 Г 


ў=1 


and the principal polarization reads 


Ре. (2а) 
EE 3 | г) 
i=l 1= 


Inspection of the above expressions shows that the values of ge (Po) and о, (Pj) 
depend on the ratio of two polynomials 5, viz. 


g — i =1 | (3) 
for Eq. (1) and (2), and 


(За) 


for Eq. (la) апа (2a). 
By means of (3) and (3a) the expressions (1) and (2) as well as (1a) and (2a) assume 
a very simple form 


25 — 1 
3—5 
еа: ©) 


The subscripts of о and Р are dropped in (4) and (5), these expressions being 
identical for fundamental and principal depolarization (polarization). 

We should like to call S the , characteristic“ of the couple of oscillators (absorption 
oscillator and emission oscillator) involved in the process of photoluminescence. 
5 is the only quantity which can be determined directly from the value of the funda- 
mental (or principal) depolarization (polarization): 

3—P 26 +1 
E TR Dee (6) 


In the case of principal polarization (the same oscillator responsible for absorp- 
tion and emission) 1 < S < 3 (cf. Eq.(3a)). In particular S = 1 for the linear oscillator, 


$ 
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ted 


S = 2 for the flat symmetrical oscillator (Г, = Г, #0 and Г, = 0), and S = 3 
for the symmetrical (isotropic) spatial oscillator (Г, = Го 20): 
In the general case (different absorption and emission oscillators), it follows from 


(3) that 1 5 S < + ә. 

For a couple of linear oscillators making an angle « 

5 = sec? a. ' (7) 

It follows from (7) that 5 = + e» (Р, = — 1/3) for х = a/2, and S = 3 for 
© == 54°44’ (Py = 0, i.e. equal to the value, which РА assumes іп the case of the 
isotropic spatial oscillator). 

The characteristic S enables the calculations of the fundamental апа principal 
depolarization (polarization) to be made very quickly. 
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ON THE INFLUENCE OF THE SHAPE OF THE NICKEL POINT 
ON THE ELECTRON IMAGE IN THE MÜLLER MICROSCOPE 


By В. Mecrewskı AND Ё. WojpA 
Department of Experimental Physics, Bolestaw Bierut University, Wroclaw 
(Received July 7, 1955) 


Previous investigations of cold emission in the Müller microscope seemed to indicate 
that the irregular distribution of the current density is mainly due to the differences of the 
work functions in different crystallographic directions and to adsorption (Müller 
1937, Becker 1951, Birkenschenkel, Haefer, and Mezger 1953, Gomer 1953). Besides, 
some authors, like Benjamins and Jankins (1940), Andreev (1952), Sokol’skaya (1952), 
Asadullin and Shuppe (1954), Spivak and Gelbert (1954) drew attention to the fact 
that changes in the electron images are caused by deformations of the point, connected 
with the migration of the atoms of the metal. The interpretation of the electron images in 
the Miiller microscope is difficult because of the fact that the above three factors appear 
simultaneously. In the present paper we intend to show that despite the accepted view 
according to which the electron image reflects the difference of the work functions, 
the shape of the point plays in certain cases a decisive and even an exclusive part. 
To achieve this aim our investigations were carried out in conditions in which the 
difference of the work functions and even the adsorption were of no importance and 
the distribution of the current density was only conditioned by the geometry of the 
point. To this end a nickel point was used, as nickel has relatively low melting point 
and the surface migration of its atoms takes place already at a temperature of 400°K. 
(Birkenschenkel, Haefer and Mezger 1953). In order to increase the migration pro- 
moting the deformation of the point observations were made on a heated point with 
a simultaneous applying of a field. While the temperature was being gradually 
raised, a decrease of the luminous areas on the screen with a simultaneous increase 
of their luminosity could be observed. This is shown in Figs. la and 1b. Comparison 
of the photograph shown in Fig. 1b with the orthogonal projection of a crystal of the 
regular system, composed of the cube faces (1 0 0), octahedron faces (1 1 1) and 
dedacahedron faces (1 1 0) (Fig. 2) permits to ascribe the dark areas to suitable pricipal 
faces and the bright önes to corners. It is easy to see that, for instance, the face (1 0 0) 
in the middle of Fig. 1b is surrounded by eight spots corresponding to the corners. 
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On the basis of our own observations as well as of those referred to in the works 
of Benjamin and Jankins (1940), Asadullin and Shuppe (1954) and Spivak апа Gelbert 
(1954) one may suppose that the corners of a nickel monocrystal heated in a strong 
electric field undergo a sharpening process. This has been confirmed by visual observa- 


Fig. 2 


tions under the electron microscope made by Gm and by Mgr Eng. A. Feltynowski, | 
as well as in the photographs given in the work of Asadullin and Shuppe (1954). 


Further investigations and an analysis of the causes of the sharpening process 
are under way. | ! 
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